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Abstract 

Given n point masses turning in a plane at a constant speed, this paper 
deals with the global bifurcation of periodic solutions for the masses, in 
that plane and in space. As a special case, one has a complete study of n 
identical masses on a regular polygon and a central mass. The symmetries 
of the problem are used in order to find the irreducible representations, 
the linearization, and with the help of the orthogonal degree theory, all 
the symmetries of the bifurcating branches. 

Keywords: ring configuration, (n + l)-body problem, global bifurca- 
tion, equivariant degree theory. 
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1 Introduction 

Consider n point masses turning at a constant speed in a plane around some 
fixed point. A relative equilibrium of this configuration is a stationary solution 
of the equations in the rotating coordinates. In this paper we analyze the 
bifurcation of periodic solutions from a general relative equilibrium. As a special 
case, we give a complete study of the polygonal relative equilibrium where there 
are n identical masses arranged on a regular polygon and a central mass. This 
model was posed by Maxwell in order to explain the stability of Saturn's rings. 
For this polygonal equilibrium, we give a full analysis for the bifurcation of 
planar and spatial periodic solutions. We will prove that, according to the value 
of the central mass, there are up to 2n branches of planar periodic solutions, with 
different symmetries, and up to n additional branches which are, if some non 
resonance condition is met, in the space with non trivial vertical components. 
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More precisely, in rotating coordinates and after a scaling in time, we look 
for 27r-periodic solutions Xj{t), for j = 0,1,..., n, where j = corresponds to 
the central body of mass while each element of the ring has mass 1, and 
frequency v. We shall prove (see the exact theorems and graphs in the core 
of the paper) that, for each integer k = 1, n, there are one or two intervals 
(one unbounded) of values of the central mass with global branches of periodic 
solutions bifurcating from the relative equilibrium. In fact, for fc = 2, ...,n — 2 
and for k = l,n — 1 and n > 6, there is an unbounded interval of values of ^ 
with two branches of planar periodic solutions (short and long period) with the 
symmetries described below. For k — l,n — 1 and n = 3,4,5,6 there is only 
one bifurcating branch in the unbounded interval, while for n = 2 there is one 
branch for any value of fi. 

For k = 2, — 2, there is a bounded interval where one has one bifurcating 
global branch for ^ < fj,k, where /Ltfc is the value of the central mass where one 
has the global branch of relative equilibria obtained in [T^. For fc = 2, ...,n/2, 
there are two bifurcating branches for < /i < toq. 

For k = n, and any n > 1, there is a global branch of planar solutions 
starting &i v = 1. 

Finally, for k = and any /i, there is a global branch of solutions 

with a non trivial vertical coordinate, which is odd in time while the horizontal 
coordinates are 7r-periodic, with the possible exception of a finite number of 
resonances at a finite number of values of /x's. 

With respect to the symmetries, the solutions in the rotating system, with 
the horizontal coordinates Uj{t) and vertical coordinate Zj{t), satisfy the follow- 
ing properties, with C = 27r/n: Mo(t) = 0, if fc and n have some common factor, 
and if not, uo(t + C) = e^*'^ ^uo(t). where k' is such that fcfc' = 1, modulo n. 

Furthemore, for j — 1, n: 

Uj{t) — e^^'^Un{t + kjC) and Zj{t) = Zn{t + kjQ). 

Except for a possible finite number of /i's and i^'s, for /x positive, bounded and 
different from ^fc, due to resonances, these solutions are geometrically different. 

By global branch, we mean that there is a continuum of solutions starting 
at the ring configuration and a specific value of the frequency, which goes to 
infinity in the norm of the solution or in the period \/v, or goes to collision, 
or otherwise goes to other relative equilibria in such a way that the sum of the 
jumps in the orthogonal degrees is zero. 

There is a vast literature on the n-body problem, many on the bifurcation 
of relative equilibria and on their stability, fewer on periodic solutions. We shall 
mention only the papers which are closer to our research, that is [26] , [32] , [27] , 
[55] . among others. The case of a satellite attracted by an array of n masses 
was treated in [IB], and previously, for the case of n = 2, by A. Maciejewski and 
S. Rybicki in ^3], using the orthogonal degree for the action of 50(2). See also 
[T7] for details. 

For our particular setting, the linearization of the system at a critical point 
is a 2n X 2n matrix, which is non invertible, due to the rotational symmetry. 
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These facts imply that the study of the spectrum of the hnearization is not an 
easy task and that the classical bifurcation results for periodic solutions may 
not be applied directly. However, we shall use the change of variables proved 
in our previous paper, [TH], in order to give not only this spectrum but also the 
consequences for the symmetries of the solutions. 

The present paper is the last part of our application of the orthogonal degree 
to similar dynamical problems. In [19], we had a complete study of the station- 
ary problem, while in [50] we have considered the analitically simpler case of 
vortices and quasi- paralel filaments. Since this paper is a continuation of jl9| . 
we shall use the results in that paper, but we shall recall all the important 
notions. 

The next section is devoted to the mathematical setting of the problem, with 
the symmetries involved in the general case and in the particular situation of the 
regular polygon. Then, we give, in the following three sections, the preliminary 
results needed in order to apply the orthogonal degree theory developed in 

, that is the global Liapunov-Schmidt reduction, the study of the irreducible 
representations, for the general and the polygonal situations, with the change of 
variables of [19], and the symmetries associated to these representations. In the 
next section, we prove our bifurcation results and, in the following section, we 
give the analysis of the spectrum in the general and the polygonal situations, 
with the complete results on the type of solutions which bifurcate from the 
relative equilibrium, in the plane and also in space. In the final section, we give 
some remarks on the resonance condition, on the stability and on the case of 
charged particles. 

Remark 1 Due to the rotational symmetry, the trivial equilibrium of the ring 
comes with an 5*0(2)- orbit and the problem is degenerate, with at least a one- 
dimensional kernel. In many previous papers with degeneracies, the authors 
look for solutions in fixed-point subspaces (odd maps, fixing the phase as in the 
classical Hopf bifurcation), in order to be able to apply simpler methods. 

In rotating coordinates, the full group is a semidirect product of with Z2 
for the general problem, and the additional group of permutations S„, for the 
polygonal relative equilibrium. 

The action ofT^ is given by time shifts and rotations in the plane, the action 
of the element ki £ Z2 is given by KiXj(t) = diag{l,—l,l)xj{—t). Thus, the 
group contains two different subgroups isomorphic to 0{2). Moreover, if one 
fixes the bifurcation map by the action of Ki , then the equilibrium must be a 
collinear central configuration, which is not the case of the polygonal equilibrium, 
or of the more general setting. 

However, the polygonal equilibrium is fixed by the action 

KXj{t) — diag{l, ~l,l)xn-j{—t), 

which is a coupling of ki with the permutation of the bodies in Dn. If one 
restricts the problem to the fixed-point subspace ofk, then one proves bifurcation 
of periodic solutions only for the blocks given by k — n and k — n/2. 
The solutions obtained for the restricted map are the following: 
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• For k = n and planar solutions, the bifurcation branch consists of solutions 
with uo{t) = and Uj{t) = e^^''Un{t), that is all the elements move as they 
were in the ring configuration. This branch was constructed in an explicit 
way in \28\/ . by reducing the problem to a (!>- dimensional dynamical system 
and a normal form argument. 

• For k = n and spatial solutions, the branch is made of solutions where the 
ring moves as a whole and the central body makes the contrary movement 
in order to stabilize the forces, that is Zj (t) — Zn (t) and zq = — ^ Zj . This 
solution is called an oscillating ring in I28f . 

• Fork = n/2 and planar solutions, the bifurcation has the central body fixed 
at the center, and two polygons of n/2 bodies pulsing each one as a whole. 

• For k — n/2 and spatial solutions, one gets the well known Hip-Hop orbits. 
This kind of solutions appears first in the paper without the central 
body. Later on, in IS81 for a big central body in order to explain the 
pulsation of the Saturn ring, where they are called kink solutions. Finally, 
there is a proof in when there is no central body. 

For the other k 's, the linearization on the fixed-point subspace of k is a com- 
plex matrix with non-negative determinant as a real matrix. This fact implies 
that one may not use a classical degree argument or other simple analytical 
proofs. See the remarks in the 4th section for a more detailed discussion. 

Note that if one wishes to restrict the problem to one of the isotropy sub- 
spaces described in this paper, one has to use anyway the orthogonal degree. Our 
approach enables us to treat the complete problem in one single system. 

Note also that one cannot ask for an arrangement of the bodies such that the 
central body remains on the x-axis and the bodies in the ring are symmetrical 
with respect to that axis, for all times, since the dynamical equations are not 
satisfied; this is why k has to reverse time. 

Remark 2 As pointed out above, the n-body problem has been the object of 
many papers, with different techniques and different purposes. For few bodies 
one may use a normal form approach, which gives good local information but 
is difficult to apply to large systems as the one we have here. The same thing 
applies to the Liapunov center theorem or the decomposition with a central va- 
riety. For the case of relative equilibria, one may consult fWj and f?j, among 
others. 

Besides the orthogonal degree, which is designed for equivariant problems 
with large symmetries, one may try to use other topological equivariant tools 
such as the equivariant Conley index or equivariant variational indices which 
effectively give a Weinstein-Moser theorem, [31, but from a finite orbit to a 
T^-orbit, which is not the case here. 

Finally, variational techniques have been quite successful in treating the prob- 
lem of existence of special solutions such as the choreographies and the hip-hop 
solutions. In particular, flO'f . UHf and classify all the possible groups which 
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give solutions which are minimizers of the action without collisions, in the plane 
for the first paper and in space for the other two. Thus, the issue is different 
from ours, since one has the proof of the existence of a solution in the large, with 
a specific symmetry, hut no multiplicity results or localization of the solutions. 

For the choreographies, following the seminal paper with no central mass, 
one has studies with more that 3 bodies and minima of the action in |^ and 
for instance. 

In the case of hip-hop solutions, the first papers such as il2j and i28j . as 
well as J^, did not use variational methods. On the other hand these methods 
were successful in 110] and 1 34 1. 

One would like to relate all the different solutions obtained by all these meth- 
ods. In the case of the analytical local solutions the linearization is also used 
in the application of any degree theory, although in the later case one may treat 
degenerate problems or large kernels. The relation between topological solutions 
and variational solutions is not easy to establish. For instance, it is one of the 
goal of JFI to prove that the bifurcation branch of hip-hop solutions may connect 
to variational solutions in the large. In particular, they give a study of the verti- 
cal bifurcation for the cases k = n and k = n/2, with no central mass and some 
conditions of non-resonance. The fact that there is no proof of the uniqueness of 
the minimizers implies that there is no rigourous proof of a global continuation 
and a connection of the hip-hop solutions to the eight choreography. 

One of the main advantages of the orthogonal degree is that it applies to 
problems which are not necessarily variational, but present conserved quantities. 
Furthermore, deformations are easy to construct. From the theoretical point of 
view the theory has to be extended to the action of non-abelian groups and to 
abstract infinite dimensional spaces ( in the case of hamiltonian systems of first 
order the degrees don't stabilize completely: see 124^ , Lemma 3.6, p. 264, ^.^d 
one has to use a simple reduction to a finite number of modes, but without 
loss of information: see the section on the Liapunov- Schmidt reduction) . See 
also ]31f for the case of gradients. Furthermore, the orthogonal degree gives 
global branches of solutions, it takes into account all symmetries and proves 
the existence of many solutions, even in the case of possible resonances, where 
the issue is to know if the obtained solution is not a subharmonic (see the last 
section or \29ij ). As pointed out above, these solutions are different. 



2 Setting the problem 

Let qj{t) G be the position of the body j <E {1, n} with mass rrij. Define 
the 3x3 matrices / = diag{l, 1, 0) and J = diag{J, 0), where J is the standard 
symplectic matrix in M^. In rotating coordinates qj{t) — e^*"^ Uj{t), Newton 
equations of the n bodies are 



J — — . — ^ ^ Uj Ui 

rUjUj + 2mj\/(jjJUj — uirrijluj — y, 'm-i'm-j — 



j=i(i#j) 
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Define the vector u as (wi, w„), the matrix of masses M as diagimil, ninl) 
and the matrix as diag{J, J). Then, Newton equations of the n bodies in 
vector form are 

Mil + 2y^MJu = W{u) with (1) 

n 

V{u) ^ — \\luj\\ + '^mimj<pa{\\uj - Ui\\). 

J = l i<3 

The function (f>a{x) is defined such that ^^(x) = — The gravitational force 

is the particular case 02 (a^) = ^/x. 

Critical points of the potential V correspond to relative equilibria of the 
n-body problem. This was part of the study done in [19] , Proposition 1. 

Remark 3 // one replaces, in the change of coordinates, the term e^*"^ with a 
complex factor ip{t) (taking qj and uj as complex functions instead of a planar 
vector), where ip satifies the equation 

ip" = -ULp/\ipf, 

the equations become 



{Mil + 2{>p)' / {ip)Mu)\'pf' = VV{u). 

In particular, the stationary solutions of this system are the same solutions 
studied in JJEl, o,nd the solutions for are those of the Kepler problem, where 
UJ is the central mass. Thus, the solutions for qj are now ellipses, parabolas 
or hyperbolas, instead of circular orbits. One may have also total collapse or 
solutions like ct^l^ , with \c\ = (^LiJ ji)^!'^ . The bifurcation analysis for periodic 
solutions starting near the relative equilibria may be performed in this slightly 
more general setting. 

To prove bifurcation of periodic solutions near a relative equilibrium, we need 
to change variables as x{t) — u{t/v). In this way, the 27r/j/-periodic solutions 
of the differential equation are 27r-periodic solutions of the bifurcation operator 

/ : i7i(R3«\vl;) X M+ ^ Ll, (2) 

f{x, v) = -u'^Mi - 2^j^uJMx + \IV{x). 

The set ^' = {a; € R'^" : Xi = Xj} is the collision set, when two or more of the 
bodies collide, and iJ|^(K^"\^') is the open subset, consisting of the collision- 
free periodic (and continuous) functions, of the Sobolev space _ff^(M^"), 

hUw'-X^) ^{x& i?L(R3") : xM ^ ^At)}- 

Definition 4 We define the action o/F = Z2 x 5*0(2) in R^" as 

p{k) ^ n and p{e) = e~^\ 

where R =^ diag{l, 1, —1) and TZ = diag{R, ...,-R), where n and 9 are elements 
ofT. 
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The group Z2 acts by reflection on the z-axis and the group SO (2) acts by 
rotation in the {x, 2/)-plane. Clearly, the potential V is invariant by the action of 
the group T. Consequently, the gradient \/V is F-orthogonal. This means that 
W is F-equi variant, p{-f)\7V{x) = \7V{p{j)x), and that \/V{x) is orthogonal 
to the infinitesimal generator, 

Aix = ■^\e=oe~-^''x = ~Jx. 

These facts can be proved directly from the definitions. 

Given that is F-equi variant and the equation is autonomous, then the 
map / is F X 5 ^-equi variant, where the action of is by time translation. 
Moreover, the infinitesimal generator of the action of S*^ in time is Aqx = x. 
Hence, the map / is F x S'^-orthogonal because of the equalities 

- \\M'/^xf - 2vV^ l^JM^'\,M"^x)^^ + V{x)\l^ = 0, 

{f{x),Jx)^, =v^ {m^/^xJm''^x) -v^ M^'^x^ + [ {VV,Jx)=Q. 

Remark 5 It is well known that the Newton equations, in fixed coordinates, are 
invariant under the action of the group of symmetries M.^ x 0(3) which gives 
the conservation laws of linear and angular momenta. In rotating coordinates, 
the potential V is invariant under the action of the subgroup Z2 x M x 0{2). 

Since the matrices J and TZ anticommute, TZS = —JIZ, then the operator 
f{x) is equivariant under the action of the full group, described by 

(T^ U KiT^) X (MU kM). (3) 

The actions of {6, ip) e and ki are given by [9, Lp)x — e~'-^^x{t + ip) 
and Kix{t) — TZix{~t), where TZi — diag{Ri, Ri) with Ri = diag{l,— 1,1). 
The action 0/ M U kR is given by the spatial reflection k, which was already 
defined, and the action ofr€Risa translation on the z-axis, rx — x + re with 
63 — (0, 0, 1) and 

e = (e3,...,e3). 

The action of the group R gives a conservation law of the linear momentum 
with respect to the spatial axis. Consequently, any spatial translation of a 
relative equilibrium is also an equilibrium. Since the orthogonal degree is defined 
only for compact abelian groups, [53], p. 70, then we have to restrict the full 
group of symmetries ([3]) to the abelian subgroup F x 5^ — T"^ x {k) . However, in 
order to follow with this procedure we have to deal with the missing symmetries 
of R by hand. The easy way to do it consists in restricting the map / to the 
orthogonal space to e G R'^". 

Definition 6 We define the map f : yVniJ|^ — > W as the restriction of the 
bifurcation map f{x) to the space 

W = {x e Ll^{R^"\-i) : {x-e)dt^O}. (4) 

Jo 
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From the previous remark, the potential V is invariant under the action of 
K. In our approach, this means that WV{x) is orthogonal to the infinitesimal 
generator of the action \/V{x) • e = 0, a fact which may be proved directly from 
the definitions. Consequently, the map / : >Vnif|^ ^ W is well defined since 

(fix) ■ e)dt = / {WV{x) ■ e)dt = 0. 



2.1 General relative equilibria 

All relative equilibria are planar. Moreover, any spatial translation of a relative 
equilibrium is also an equilibrium. Nevertheless, the only relative equilibrium 
in W of the family of translations is the one in the {x, t/)-plane. Therefore, the 
positions {uj , 0) form a relative equilibrium in W if they satisfy the relations 

n 

rrii , -I with U j G K . 

Remark 7 Since the potential is homogenous, then any scaling of a relative 
equilibrium is also an equilibrium. Hence, in principle we may choose the relative 
equilibrium with cj = 1. However, we shall leave u as a parameter because it is 
easier to fix the norm of the polygonal equilibrium than the frequency. 

Since every relative equilibrium xq is planar, then the action of k e F leaves 
fixed Xq. Therefore, the isotropy subgroup of in F x is Fa;^ x with 
= (^) = ^2- Thus, the orbit of xq is isomorphic to the group S^: the orbit 
of the equilibrium xq consists of all the rotations of xq in the (x, ?/)-plane. 

Notice that the generator of the spatial group F at xq is AiXq ~ ~Jxq, then 
—Jxq is tangent to the orbit and must be in the kernel of D'^V{xq). 

Now, since Z2 is in the isotropy subgroup of xq, then D'^V{xq) is Z2- 
equivariant. From Schur's lemma, see [24j . p. 18, the matrix D'^V{xq) must 
have a block diagonal form. In the following proposition we prove directly this 
fact. 



Proposition 8 Let Aij be the 3x3 blocks such that 

D^Vixo) = {A,,)7j=,, 
then the matrices Aij have the diagonal form 

Aij = diag{Aij,a.ij). 

Let dij be the distance between Ui and Uj, and let (xj,yj) be the components 
of Uj. For i 7^ j, we have that aij = mimj/d"^^ and 

m,mj ( {a + l){x,~Xjf -d}^ {a + l){x^ - Xj){y, - y,) \ 
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Moreover, for i — j we have the equalities 

n 

a-ii = - ^ aij and Ai^ = Lorriil - ^ Aij. 

Proof. The function (dij ) has the second derivative matrix 

a + 1 [ (xi-Xj)^ (x, - Xj){y, - yj) 

Dl,(l}a{drj) =^ \ {x,-Xj){y^-Vj) iVt-yj)^ , . 

% \ / % 

For i ^ j, since V„,0(dij) = -V„^0((iij ), then 

Aij = mimjDujV u,(t>{.dij) = -mimjDl^(p{dij) = diag{A,j , aij). 

And for i = j the matrix satisfies 

Ail = ujmj +^^mimjDl^(f>{dij) = ujmj -^^Aij = diag^Au^au). 



If we analyze the n-body problem in the plane instead of the space, then the 
potential of the planar problem has the Hessian D'^V{xo) = {Aij)2j^i. Because 
of this fact, we say that {Aij)^^^^ will give the planar spectrum and {aij)2j^i 
the spatial spectrum. 



2.2 The polygonal equilibrium 

Identifying the real and complex planes, a relative equilibrium is given by the 
positions {aj,0) G K'^ with Uj G C. Defining ( = 27r/n, the polygonal equilib- 
rium is formed by n + 1 bodies: one body, with mass mp — /i, at oq — 0, and 
one body for each j g {1, ...,n}, with mass rrij = 1, at aj = e'-'''. 

It is easy to prove, for instance jl9| . Proposition 1, that the ring configuration 
d = (ag, a„) is a relative equilibrium when lu = fi + si, with 



In this paper, the norm of the polygonal equilibrium is fixed and oj = fj.+si is 
a free parameter. In the paper [TH], Theorem 24, we have proved bifurcation of 
relative equilibria from the polygonal equilibrium with the parameter fi . Now, 
we wish to analyze the bifurcation of planar and spatial periodic solutions. 

Definition 9 Let Sn be the group of permutations o/{l,...,n}. We define the 
action of Sn in R2("+i) as 

p{j){xo,xi,...,Xn) = (xo,a;-y(i), 
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As for the general case, the map / is (Z2 x SO{2)) x 5" -^-orthogonal. Moreover, 
since n of the bodies have equal mass, then the gradient W is S'„-equivariant. 
Therefore, we may consider the map / as F x S'^-orthogonal with respect to the 
abelian group 

r = Z2 X Z„ X S0{2), 

where Z„ is the subgroup of S'„ generated by ({j) = j + 1. 

Let Z„ be the subgroup of F generated by (C, C) G x S0{2) with ^ = 
27r/n € 50(2). Since the actions of (C, C)i given by taking the j-th element in 
the plane xj to e~''''Xj+i, for j e 1, ., ., ., n and Xq to e~''''Xo, and k G Z2, fixing 
the plane, leave the equilibrium a fixed, then the isotropy group of a is Fg x 
with 

Fa = Z2 X Z„. 

3 The Liapunov- Schmidt reduction 

The orthogonal degree was defined only in finite dimension. There is a the- 
oretical difficulty for the extension to abstract infinite dimensional spaces, in 
particular, for strongly indefinite problems, in the sense that the degrees de- 
fined on finite dimensional approximations do not stabilize, as the classical 
Leray-Schauder degree, when increasing the dimension, although the changes 
are easy to compute: see, for instance, the case of hamiltonian systems treated 
in |24j . Lemma 3.6, p. 264. In the present situation, one has a very simple way 
of overcoming this difficulty, without any loss of information and, at the same 
time, avoiding many technical problems of functional analysis. 

Thus, we must perform first a reduction to some finite space . 

The bifurcation operator /, in Fourier series, is 

fix) = Y^il-'v^Mxi - 2lu^{iJ)Mxi + gi)e''\ 

where xi and gi are the Fourier modes of x and W{x), which is relatively 
compact, by Sobolev embedding, with respect to the second order terms. Since 
all the masses are non-zero, then the matrix l^i/^Ail — 2ilv^fbjJ M. is invertible 
for all large /i^'s, for v > Therefore, if we choose a big enough p, we can 
solve the modes |/| > p in terms of the remaining 2p + 1 Fourier modes, on any 
bounded subset of W with frequencies v uniformly bounded from below, with 
an application of the global implicit function theorem, [24], p. 257. 

In this way, we get that the zeros of the bifurcation operator / correspond 
to the zeros of the bifurcation function 

f{xi,X2{xi,i^),iy) = Y.ifiy^Mxi-2liyV^{iJ)Mxi+gi)e'^\ 

where Xi corresponds to the 2p+l Fourier modes and X2 to the solution for the 
remaining modes. Then, the linearization of the bifurcation function at some 
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equilibrium xq is 

\i\<v 

Since the bifurcation operator is real, the hnearization is determined by the 
blocks M[lv) for I S {0, where M(y) is the matrix 

M{v) = v^M - 2v^{iJ)M + B^Vix^;). (6) 

Each one of the blocks Milv) corresponds to a Fourier mode:since the orbits 
are real, one has that the a;;, for negative /'s are the conjugates of X-i. This is 
why we are considering xi € C, for / € {0, Since we restrict the operator 

](x) to the subspace of iJ|^(M^"\^) orthogonal to e, then the block for the first 
Fourier mode Af(0) must be restricted to the orthogonal space to e. We shall 
denote this restriction as M(^')^^ . 

4 Irreducible representations 

The following step in order to use the orthogonal degree, [53] Theorem 3.1, p. 
247, consists in identifying the irreducible representations, that is, for the action 
of the isotropy group Fg and obtaining the decomposition of the matrix M(v). 

4.1 The general relative equilibria 

We show first, for n arbitrary bodies in the plane, how the matrix M (y) decom- 
poses into two blocks by the action of the isotropy group Fa;^ = Z2. The two 
blocks correspond to the planar and spatial spectra as in the previous decom- 
position of D'^V{xq). 

Since the group Z2 acts by reflection on the z-axis, then the equivalent 
irreducible representations for C'^", leading to the decomposition of M(v)^ are 

Vb = {(2;i,2/i,0,...,x„,y„,0) : Xj.y.j £ C} and 
= {(0,0,zi,...,0,0,z„) eC}. 

The action of k e Z2 on Vq is given by p{k) = I and on Vi by p(k) = — /. Define 
the isomorphisms To : C^" Vq and Ti : C" Vi as 

To{xi,yi, ...,Xn,yn) = (a;i,yi,0, a:„, 0), 
Ti{zi, ...,z„) = (0,0, zi, ...,0, 0,z„), 

and define the linear orthogonal map P 

P{xi,yi,zi,...,Xn,Vn, Zn) = To(a;i,yi, ...,x„,?/„) -I- Ti{zi, ...,z„). 

The transformation P rearranges the planar and spatial coordinates. 
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Since Vq and Vi arc subspaccs of equivalent irreducible representations, 
Schur's lemma implies that the matrix M{i') must satisfy 

p-^M{iy)P = diag{Mo{ty),Mi{iy)), 

where Mq and Mi arc respectively 2n x 2n and n x n matrices. In fact, we will 
exhibit explicitly the matrices Mq and Mi, in the following result: 

Proposition 10 Define the matrices M.i = (mi, m„), A^2 = diag{mi,mi, ...,mn,mn), 
and J = diag{J, J), then the blocks Mq and Mi are 

Mo{v) = v^M2 - 2v^{iJ)M^ + {A^)l^i, 
Mi{v) = v'Mi + (a,,)%^^. 

Proof. Since P rearranges the planar and spatial coordinates and D'^V{xq) = 
(Aj)rj=i. with Aij = diag{Aij,aij), then 

p-'D^V{xo)P = dzag{{A,j)-^„ (a,,)^;=i). 

Moreover, 

p-^MP = diag{M2,Mi) and P'^JP = diag{J,0). 

Therefore, 

p-^M{u)P = u^diag{M2,Mi) - 2vs/Zjdiag{{iJ)M-i^,Q) 
+ diag{{Aij)l^i,{aij)l^i). 

■ 

The block Af (0) must be restricted to the space orthogonal to e. Since 
e G Fi, with Ti(l,...,l) ~ e, then, in the new coordinates, the block Mi(0) 
must be restricted to the space orthogonal to (1, 1). 

For the study of the planar n-body problem, this procedure gives the matrix 
Mo{iy) instead of M{iy). 

PlanEir representation 

The action of the isotropy group Z2 x on the space Vq is given by 

(k, ip)x = e^'^x. 

Since k G Z2 fixes the points of Vq, then the space Vq has its isotropy subgroup 
generated by k, 

Z2 = (k) . 

Thus, solutions x{t) with isotropy group Z2 must satisfy 

x{t) = Kx{t) = Rx{t). 
Consequently, orbits in Vq are planar, Zj{t) = 0. 
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4.1.1 Spatial representation 

In Vi the action of the group Z2 x is 

(k, ip)x = —e^'^x. 

Since the action of (k, tt) fixes the points of Vi , then the space Vi has its isotropy 
subgroup generated by (k^tt), 

Z2 = ((k, tt)) . 

Solutions x(t) with isotropy group Z2 must satisfy 

x{t) — (k, T:)x{t) — Rx{t + tt). 

Consequently, orbits in Vi satisfy the symmetry 

(t) = (t + T^), Vj [t] = V] {t + tt) and zj (t) = -zj (t + tt) . (7) 

From the symmetry (O, we have that Zj[tj) = for some tj. This means 
that these solutions oscillate around the (x,y)-plane and that the projection 
of the curve in the (x, ?/)-plane is 7r-periodic. Furthermore, these solutions go 
around the projected 7r-periodic curve once with the spatial coordinate Zj{t) 
and once with —Zj{t). These solutions look like eights around the equilibrium 
points, and we shall call them spatial eights. 

4.2 The polygonal equilibrium 

In the case of the polygonal equilibrium, the isotropy group has also the action 
of Z„. Notice that the action of (C, C) G in Vfc, fc = 0, 1, is given by 

(C,C)7o(a;o,2/o, ...,2:„,2/„) = To(e-^'^(xo, yo, ••■,yc(n)))' 

(C, 0^1(2:0, Zn) = Ti(zo, 2;c(n)). 

Thus, we expect a decomposition of the space Vk into smaller irreducible repre- 
sentations due to the action of Z„ . 

4.2.1 Planar representations 

We give first the decomposition of the space Vg into smaller irreducible represen- 
tations. Since the representations for n = 2 are different from those for n > 3, 
we shall restrict the study to the case n > 3 and we shall give some comments 
on the case n = 2 at the end of the paper. 

Definition 11 For k £ {2, ...,n — 2,n}, define the isomorphisms : — >■ Wk 
as 

Tk{w) = {0,n-^/^e^''''+'^^w,...,n-^/^e"^''''+'^^w) with 
Wk = {(0, e^^'^^+'^^'^w, e"(^'=^+'^)'^w;) : w £ C^}. 
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For k £ — 1}, define the isomorphism Tk : ^ Wk as 

Tk{a,w) = (vfca,n-i/2e(^'=-f+^)««;,...,n-i/2e"(''^-f+-^)'^u)) with 
Wk = {{vka, e^^'^^+'^^'^w, e"(*'=^+'')^w;) ■.aeC,w£ C^}, 
where vi and are the vectors 

vi = 2^1/2 (^i^ ^ 2-^/2 (1, -i) . 

We have proved, in the paper [TO], p. 3207, that the subspaces Wk form the 
irreducible representations of Z„ in C^*^"+^-' ~ Vq. Also, we proved, in [T5] . 
Proposition 5, that the action of (C, C, (^s) G Z„ x S*^ on Wki, for the /-th Fourier 
mode, is given by 

Consequently, the isotropy subgroup of Fg x §^ for the block Wk , for the fun- 
damental mode, is generated by k G Z2 and {(, (, — fcC) € x S*^, that is, 

Z„(A:)xZ2 = ((C,C,-fcC)>x(K). 
Since the subspaces Wk are orthogonal to each other, then the hnear map 

Pw = y^rfc(wfc) 

fe=l 

is orthogonal, where w = {wi, w„), with e C"^ for k — 1. n— 1 and Wk G 
for the remaining /c's. 

Since the map P rearranges the irreducible representations of Z„, from 
Schur's lemma, we must have that 

P-\Aj)7j^,P = diag{Bi,...,B„), 

where Bk are matrices such that {Aij)fj^^Tk{w) — Tk{Bkw). In fact, in the 
paper [Tn|, Propositions 18 and 19, we gave the blocks Bk for the general case 
a > 1. In the next proposition, we state this result. 

Proposition 12 Set a+ = {a + l)/2 and a_ = (a — l)/2. Define 

Uk = a^isk+i + Sk-i)/2, Pk = a+{sk ~ si) and jk = a_(sfc+i - Sk-i)/2 

with 

_ 1 sin^(fcjC/2) 
2«^^sin«+i(iC/2)' 

Then, the blocks satisfy Bn-k = Bk and they are given by 

Bk =a+{I + R)^i + (si + ak)I - PkR - IkU , 

for k £ {2, n — 2, n}, and 

II (si + + na^) — ■\/n/2/ia — •y/n/2/ii 

Bi = I —^Jripiiiot si + ai + (a + a\i 
\J n/2^,i —aii si + a\ 
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Since Mi{u) is Z„-equivariant, from Schur's lemma, we have that 

P~'^Mo{u)P = diag{moi{i^), ...,mon{'^)), 

where the matrices mofc(z^) must satisfy MoTfc(w) = Tfc(mofcw). 

Proposition 13 The matrices mok{v) satisfy mofe(i^) = 'mo(n-k){~^) with 

TOofe(z^) = f"^! - 2vy/oj{iJ) + Bk 

for k e {2, n — 2, n}, and 

iTioiiv) = diag{iJL, I) — 2v^/ujdiag{n,iJ) + Bi. 

Proof. The matrix A42 is diag{iJ,, fi, 1, 1). For k £ {2, n—2, n}, the matrix 
M2 satisfies M2Tk{w) = Tk{w). Since {iJ)Tk{w) = Tk{iJw), then 

MoTkiw) = (1/^X2 - 2iyV^M2{iJ) + (^i,)"j=i) Tk{w) = Tk{mokw). 

For k = 1, the matrix A42 satisfies A42Ti{w) = Ti{diag{iJ,, l,l)w). Since 
{iJ)vi = vi, then {iJ)Ti{w) = Ti{diag{l,iJ)w). Therefore, 

MoTi{w) = {u^M2 - 2v^M2{iJ) + (^y)S=i) ^^iH = Tx{mmw). 
For fc = n — 1, using {iJ)v2 = —V2 and a similar argument, we may prove that 

mo(„_i)(i') = v^diag{^, I) — 2v-Jbjdiag{—^,,iJ) + Bn-i- 

The equalities 'tiq(„_j,-) (i^) = rhoki—i^) follow from -B„,_fc = B^- ■ 

Remark 14 In fact, the equality mo(„_fe)(i^) = Tfn^k{—v) is a consequence of 
the Za2 -equivariant property of the operator by the action KXj{t) = TZiXn-j{—t), 
for j = 1, n — 1, and KXj{t) = TZiXj{—t), for j = and j = n. 

In fact the bifurcation operator f{x) is equivariant under the action of the 
group 

(Z„ X T^) X k(Z„ X T2), 
where the action of k is given, for the n equal mass bodies, modulus n, by 

KXj{t) = R\Xn-j{—t) with R\ = diag{l, —1, 1), 

while, for the central body, j — 0, the action is 

KXo{t) = RiXo{—t). 

Let Zk be the coordinate for the block mok-The action of k on {zk,Zn-k) is 
given by 

k{Zk,Zn-k) = {RZn-k,RZk), 

where R = diag{l,—l) for k ^ {l,n — 1} and R = diag{l, 1,-1) for k G 
{1, n — 1}, and k{zn) = Rzn- 
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Now, since the polygonal equilibrium a is fixed by k, then the lineariza- 
tion of the bifurcation map f{x), represented by the blocks m^k, must be {k) 
-equivariant. This means that 

kdiag{mok,mo{n^k))izk, Zn^k) = diag(™ofe, ™o(n-fc))'«(^fe, ^n-fc)- 

From the first coordinate one gets the equality RrfiQ(^n-k) = mokR, or equiv- 
alently 

ruQk Rmo(^n-k)R, 

which is equivalent to the equality moki^') — fho{n-k){^^)- 
Now, if we restrict the bifurcation operator f{x) to the fixed-point space of 
k, then we get that the coordinates Zk of the blocks wofe are related by 

Zn-k = Rzk and z^ = i?z„. 

Using the isomorphism Tzk — {zk,Rzk) for k (f. {n/2,n} and T{xk,yk) = 
{xk,iyk) € for k G {n/2,n}, one may see that Afo(i^), in the fixed-point 
space of k is equivalent to the matrix 

(moi,mo2,TOo3, ■•■,TO0(n/2),™0n) 

defined in the space 

X X X ... X R2 X 

Applying degree in the fixed-point space of k, one proves bifurcation of pe- 
riodic solutions for the blocks given by mon and niQn/2- Moreover, since the 
action of k and 

(C,C,0) e Z„ X X 

on is trivial, as well as the action of k and (C, C, tt) on mo„/2, then one 
can deduce also that these equilibria will have the symmetries that we described 
above. 

The remaining blocks rriofc for k (f. {n/2,n} are defined in complex subspaces, 
then they have always non negative determinants, as real matrices. Thus, there 
is no change of a standard degree in the isotropy subspace. Analytical studies 
with normal forms of high order and additional hypotheses of non-resonance are 
proposed in for these cases. One could also use the gradient structure and 
apply the results for bifurcation given in \21}/ . p. 100, based on Conley index. 
However, these results do not provide the proof of the existence of a continuum, 
something which follows from the application of the orthogonal degree. 

4.2.2 Spatial representation 

We analyze now the decomposition of Vi . 



16 



Definition 15 For k e {1, ...,n — 1}, we define Tk ^ Wk as 

Wk = {(0, e''"^w, e"''=^w) :weC}. 

And, for k = n, we define r„ : — )• W„ as 

Tn{a,w) = {a,n^^^^w, ...,n^^^^w) with 
Wn = {{a, w, w) : a,w G C}. 

The action of {(, () in the subspace Wk is given by 

{C,OTk{w)=Tk{e'''<w). 

Consequently, the decomposition in irreducible representations of the space 
£n+i ^ £qj. |.]-^g group Z„, are the subspaces Wk- The actions of the el- 
ements K e Z2, (C, C) € and (fi G in Wk are 

p{^) = -l,p(C,C) = e^"^ and p(^) = e''^. 

Since the elements (CiC)~^C) G Z„ x 5^ and {k,tt) £ Z2 x act trivially on 
Wk, then the isotropy group of Wk is generated by (C, Ci "'^'C) and (/t, tt), 

Z„(fc)xZ2 = ((C,C,-fcC))x((K,7r)). (8) 

Since the spaces Wk are orthogonal to each other, then the linear map 

n 

Pw = Y,Tk{wk) 

is orthogonal, where w = (wi,...,w„), with w„ e and Wk € for the 
remaining fc's. 

Since Mi(:/) is Z„-equivariant, by Schur's lemma, the matrix Mi(i^) must 
satisfy 

p-H4i{iy)P = diag{mii{iy), mi„(i^)), 
where the blocks mifc(i/) are such that MxTkiw) = Tk{mikw) for k £ {1, ...,n}. 

Proposition 16 We have mife(z/) = v'^ — + Sk) for k e {1, n — 1}, and 

Proof. Wc denote the coordinate Wi G of the vector w = {wq, ...,Wn) by 
[w]i. For k e {1, ...,n — 1}, if / 7^ 0, from the definition 

n 
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Since aij = minij/dij, with df^ = 4sin^((Z — j)C/2) for l,j G then 
O'lj = <in{j-i): with (j — I) € {1, n}, modulus n. From the equahty aijC^^*"^ = 
6*""^ (a„(j_i)e'(-'~')'"^), we have that 

[{a,,)Tk{w)]i = n-i/2g.ifcC a„,e^^'=^ j - [Tk{bkw)]i , 

where bk is the sum between parentheses. 

In order to calculate bk, notice that a„„ — — Y^^Zo '^"i ^^'^ "^"O = then 

n-l 

Since a„j = then 

y (e^^'^ - 1) a„, = - y '^^"'^y'^ = 



-2"+isin"+^(jC/2) 



Thus, we have {aij)Tk{w) = Tk{bkw), with bk = -(^ + Sfc). 
If ? 7^ 0, then, 

n 

[{aij)Tn{a, w)]i = aioa + rr^^"^ ^ ajj-w, 

with a, w € C. Using the equalities a;o = /i and X]j'=i '^ij — ~'^io, we have 

[{aij)Tn{a,w)]i ^ n~^^^{^/nna - ^iw). (9) 

If ? = 0, then 

n 

[(fly )T„(a, i(;)]o = aooa + n^^^'^ ^ aojw;. 

Using the equality ooo = —n^i, we have 

[{aij)Tn{a,w)]o = —n^a + ^/nfjLW. (10) 
From the equalities ([9]) and JTU]), we conclude that 



Here the matrix Mi is diag{fi, 1, 1). Then, the statements of the propo- 
sition follow from the fact that A4iTk{w) ~ Tk{w) for k E {l,...,n — 1}, and 
MiTn{a,w) ^ Tn{fJ-a,w). ■ 
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Since the block Mi(0) is restricted to the space orthogonal to (1, 1), and 
since Tn{l,\/n) = (1,...,1), then the block mi„(0) must be restricted to the 
subspace orthogonal to {1,^/n). 

As we expected, the eigenvalues of mi„(0) are and — (n + with eigen- 
vectors {l,^/n) and {^/n,—l) respectively. Therefore, the matrix mi„(0), on 
the subspace orthogonal to (1, y/n), is equivalent to 

mi„(0)-L = -(7i + l)/i. 

The map rni„(0)"'" is invertible for n 0- For fi = the central body has zero 
mass, and, in this case, the coordinate corresponding to min(O)"'" may be taken 
away. 

Remark 17 In the case of the spatial blocks, let Zk be the the coordinate of the 
block niik-The action on this coordinate is given by 

for k — 1 , . . . , n — 1 and kzn — Zn- 

Since the blocks mik must be (k) -equivariant, then one gets the equality 

mik = mi(„_fc). 

Thus, when restricting the bifurcation operator f{x) to the fixed-point space 
of k, the coordinates Zk of the blocks mik o,re related by — Zk, and z„ — z„. 

Using the isomorphism Tzk = {zk,Zk) for k ^ {n/2,n}, the inclusion of W 
in C for k = n/2, and the inclusion o/M^ in for k — n, one may prove that 
Mi{iy) in the fixed-point space of k, is equivalent to the matrix 

(mil, mi2, TOl(n/2), »7ii„), 

defined in the space 

C X C X ... X M X 

Applying a degree argument in the fixed-point space of k, one may prove 
bifurcation of periodic solutions for the blocks given by mi„ and TOi„/2 . The 
remaining blocks mofc for k ^ {n/2,n} are defined on complex subspaces, then 
they have always non-negative determinants, as real matrices. Using the com- 
ments of the previous remark, one may give a different proof of the results on 
vertical bifurcation discussed in JS^. 

5 Symmetries 

Before we show the existence of bifurcation points, we wish to describe the 
symmetries of the solutions. That is, solutions with isotropy group Z„(fc) x Z2 
for the blocks mofc(z/) and Z„(fc) x Z2 for the blocks mik{v). 

The coordinates {xj,yj, Zj), of the body j G {0, n}, will denote the com- 
ponents of the solution x{t). 
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5.1 Planeir solutions 

In this part wc describe the symmetries of the isotropy group Z„(fc) x Z2. As 
we have seen, due to the group Z2, these solutions must be planar. We shall 
identify the real and complex planes by uj — xj + iyj . 

Now, since (C, Ci "'''C) € Z„(A;), the solutions with isotropy group Z„(fc) must 
satisfy 

Uj{t) = e-"^u^^j)it-kO. 

Remark 18 Notice that, ifuj{t) is a solution, forv, with symmetry Zn^k) , then 
Uj{—t) is a solution, for —v, with symmetry Zn{n — k). In fact, the bifurcation 
phenomena are related by the fact that mo(„_fe)(i') = mofe(— z^)- 

In order to describe the symmetries of the group Z„(A:) we need the following 
definition. 

Definition 19 For each fixed k, let h be the maximum common divisor of n 
and k. We define 

n = n/h and k = k/h. 

For the central body uq one has the following symmetries. 

Proposition 20 If h > 1, the central body remains at the center uo{t) = 0. // 
h = 1, the central body satisfies 

uo{t + C) = e-"''^uo{t), 

were k' is such that k'k = 1, modulus n. 

Proof. Since C(0) = 0, the central body has the symmetry 

uoit) = e-'^'^Ko{t-k{lO)- 

For h> 1 take I = n. Since fc(n(^) = 27rfc and n( = 27r/h, then the central body 
satisfies 

uoit) = e-'("f)wo(t - k{nC)) = e-''^'^^'^\o{t) = 0. 
For h = l take I = k', then the central body satisfies 

uo{t) = e-'('='^)wo(i - k{k'0) = e-'^'^'^uoit - ()■ 

■ 

In order to describe the symmetries of the n bodies with equal masses, we 
use the notation Uj = Uj+kn, for j G {1, ...,n}. Then, ({j) = j + 1, and the n 
bodies satisfy 

Uj+i{t) = e'^'^ui{t + jkO. 

Thus, each one of the n bodies follows the same planar curve, but with a different 
phase and with some rotation in the {x, y)-plane. 
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Now, let us show some examples of these symmetries. 
For Z„(n), the symmetries arc 

uo{t) = and = e *-'^-ui(t). 

Thus, the central body remains at the center, and the other n bodies form a 
n-polygon at any time: see the figure for n = 3. 



(a) Symmetries of Z„(l). (b) Symmetries of Z„ (3). 

Figure 1: For n = 3. 

For Z„(l), the symmetries are 

uo{t + C) = e-'^uo{t) and Uj+i{t) = e '^^ui{t + jQ. 

Therefore, the central body is determined by the time interval [0,^), and the 
rest of its orbit is given by rotations. The other n bodies follow the orbit of 
one of them, but with a synchronization between the change of phase and the 
(a;, t/)-rotation: see the figure for n = 3. 

For a group Z„(A;), with h = l, the symmetries are 

uo{t + C) = e-''='^uo(t) and Uj+^{t) = e '^^ux{t + j{kQ). 

These symmetries are similar to the ones for Z„(l), but now there is a permuta- 
tion between the phase and the planar rotation. For instance, one may compare 
the cases Z5(l) and ^5(2), for n = 5. 

The effect of the group ^5(1) was already described. In order to show the 

symmetries for the group ^5(2), notice that /c' = 3, then the central body visits 
the points: uq{Q = e~^^'^UQ, Uo(2C) = e~*'»UO) '"o(3C) = e~^*^uo and xo{AC,) = 
e~^'^uo. The other n bodies have orbits: e + 2^), e + 4C), 

e *3C„^(t + Q and e '4«Mi(t + 3C). 

For the group 1in{k) with h> 1, the symmetries are 

uo{t) = and Uj+i{t) = e ^^'^ui{t + jk{2-K/n)). 
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(a) Symmetries of Z„(l). (b) Symmetries of Z„(2). 



Figure 2: For n = 5. 




Figure 3: Symmetries of Z„(2) for n = 4. 

In this case the central body remains at the center. Moreover, the other n bodies 
follow the same curve with a different phase and the planar rotation determined 
by multiplying by k in Zfj. See the example n = 4 and k = 2 . 

Remark 21 In fixed coordinates, the solutions areqj{t) = e'^^*Uj{vt). Reparametriz- 

ing the time, we have that qj{t) = e'^*/''uj(f), where Uj is a 2-1: -periodic func- 
tion. Set Q. = \ — ky/uj/v, then, for j G {1, ■■■,n}, we have that 

9,■+l(^) = e"^/'^^^,+l(^) 

= e *^«e^*^/''ui(t + jkQ = e ''^^\i{t + jkQ. 

In particular, when the central body has mass zero, we are considering the 
n-body problem with equal masses. If G nZ, the solutions with isotropy group 
I'n{k) satisfy 

qj+i{t) = qi{t + jkC). 
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These are solutions where all the bodies follow the same path, and they are 
known as choreographies, J^. 

5.2 Spatial solutions 

Now we wish to describe the symmetries of a sohition with isotropy group 
Z„(fc) X Z2. As we saw for the group Z2, these solutions satisfy 

Uj{t) — Uj{t + tt) and Zj{t) — —Zj^t + tt), 

where Uj — Xj + iyj is the projection in the (a;, y)-plane. 

Since the group Z„(A;) is generated by (C, C: — ^C): these solutions satisfy the 
symmetries Uj{t) = e~*^it^(j)(t — kQ and 

zAi) = HU)ii-^0- (11) 

Therefore, the projection in the (a;, ?/)-plane follows a 7r-periodic curve with 
the symmetries of the previous section. The spatial solutions go around this 
projected curve twice, once with Zj (t) and once with Zj {—t). Hence, the solutions 
look like spatial eight where the spatial positions satisfy the symmetries (jlip . 

Hence, we need to identify the spatial symmetries pT|) . For the central body 
we have that 

Zo(t) = Zo{t + kC). 

Using the notation Zj = Zj^kn for the other n bodies, their spatial positions are 
related by 

Zj+i{t) = zi{t + jk(^). 

Therefore, the spatial curves of the n bodies are determined by just one of 
them. To see one example, we suppose that n — 2m and choose k ^ m. In this 
case the central body remains at the center. Moreover, the n bodies with equal 
masses satisfy 

Uj+i{t)^e'^'^u,it + j7r)=e'3'^u,{t) 

and 

z,+i{t) ^ ziit + jn) = i-iyz,{t). 

Thus, there are two m-polygons which oscillate vertically, one with zi{t) and 
the other with —zi{t). Furthermore, the projection of the two m-polygons in 
the plane is always a 2m-polygon. These solutions are known as Hip-Hop orbits. 

Remark 22 An element in the intersection of two isotropy subspaces, for ki 
and ^2, will be such that Uo{t) — 0, Uj{t), ZQ{t), Zj{t) will be (fc2 — ki)C,-periodic, 
that is, if h is the greatest common factor of n and k2 — ki and fih = n, then 
the period will be 27r/n. 



23 



6 The bifurcation theorem 



The orthogonal degree is defined for orthogonal maps which are non-zero on the 
boundary of some open bounded invariant set. The degree is made of integers, 
one for each orbit type, and it has all the properties of the usual Brouwer degree. 
Hence, if one of the integers is non-zero, then the map has a zero corresponding 
to the orbit type of that integer. In addition, the degree is invariant under 
orthogonal deformations that are non-zero on the boundary. The degree has 
other properties such as sum, products and suspensions, for instance, the degree 
of two pieces of the set is the sum of the degrees. The interested reader may 
consult Chapters 2 and 4,[T] and jSl, for more details on equivariant degree 
and degree for gradient maps. 

Now, if one has an isolated orbit, then its linearization at one point of the 
orbit xo has a block diagonal structure, due to Schur's lemma, [24], Lemma 
7.2, p. 30, where the isotropy subgroup of xq acts as Z„ or as . Therefore, 
the orthogonal index of the orbit is given by the signs of the determinants of 
the submatrices where the action is as Z„, for n = I and n = 2, and the 
Morse indices of the submatrices where the action is as S^, [21], Theorem 3.1, 
p. 247. In particular, for problems with a parameter, if the orthogonal index 
changes at some value of the parameter, one will have bifurcation of solutions 
with the corresponding orbit type. Here, the parameter is the frequency i^, [24], 
Proposition 3.1, p. 255. 

For a fc-diniensional orbit with a tangent space generated by k of the in- 
finitesimal generators of the action of the group, one uses a Poincare section for 
the map augmented with k Lagrange-like multipliers for the generators. (See 
the construction in [21], Section 4.3, p. 245.) For instance, for the action of 
50(2), the study of the zeros of the equivariant map F{x), orthogonal to the 
generator Ax, is equivalent to the study of the zeros of F(x) + XAx, if x is not 
fixed by the the group, i.e., if Ax is not 0, for which A is 0. In this way, one has 
added an artificial parameter. This trick has been used very often and, in the 
context of a topological degree argument, was called "orthogonal degree" by S. 
Rybicki in [30 . See also [11^ and [22], p. 481, for the case of gradients. The 
general case of the action of an abelian group was treated in [23] . The complete 
study of the orthogonal degree theory is given in f5?^. Chapters 2 and 4. 

Any Fourier mode will give rise to an orbit type (modes which are multiples 
of it), hence one has an element of the orthogonal degree for each mode. Fur- 
thermore, if x{t) is a periodic solution, with frequency then y{t) = x{nt) is 
a 27r/rt-periodic solution, with frequency v/n. Hence, any branch arising from 
the fundamental mode will be reproduced in the harmonic branch. If one wishes 
to study period-doubling, then one has to consider the branch corresponding to 
TT-periodic solutions, [24 , Corollary 2.1, p. 219. 

The complete study of the orthogonal degree theory is given in [24] , Chapters 
2 and 4. 
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6.1 General relative equilibrium 

To use successfully the results of the orthogonal degree from [21], Proposition 
3.2, p. 258, we need to veriiy the hypothesis that the orbit Tx is hyperbolic near 
a bifurcation point, [23], Definition 2.2, p. 222. In particular, we need to prove 
that the kernel of M(O)-'-'^ is generated by AiXq ~ —Jxq- Since AiXo S Vb 
and e S Vi, an equivalent condition consists in showing that the kernel of 
Mo(0) is generated by Tq^{Aixo), and that the kernel of Mi(0) is generated by 
T-i(e) = (l,...,l). 

Definition 23 Following JSJI, p. 258, we define a as the sign of the determi- 
nant of Mq{0) in the space orthogonal to Tf^^{j'xo). 

Now, if (zi, 2„) is in the kernel of Mi(0) — (aij), then Zj = Zi, because 

n n 

Consequently, the kernel of Afi(O) is generated by (1, 1), and the hypothesis 
of [23], Proposition 3.2, p. 258, near a bifurcation point, is assured by the 
condition cr 7^ 0. 

Definition 24 Let nj{i') be the Morse number of Mj{v), for j = 0,1, and 
define 

Vjiy) = ^i.'^ji^ - P) - + p)). 

The number rijiyo) represents the change of the Morse index at the point 
Uq. Applying the bifurcation theorems of ([24]), Remark 3.5, p. 259, we get the 
following result: 

Theorem 25 If rij{vj) is different from zero, the equilibrium xq has a global 
bifurcation of periodic solutions from 2T:/i>j, with isotropy group Z2, for j = 0, 
and Z2 , for j — I. 

The bifurcation branch is non-admissible when a) the norm or the period 
goes to infinity or b) the branch ends in a collision path. In any other case, 
we say that the bifurcation is admissible. By global bifurcation we mean that: 
if the branch is admissible, then the branch must return to other bifurcation 
points, and the sum of the local degrees at the bifurcation points rjk{i') must be 
zero. See [21] for this kind of arguments. 

6.2 The polygonal equilibria 

We have decomposed the matrix Af(i^) into the two blocks Mq{i^) and Afi(t^). 
Thus, there is a bifurcation of periodic solutions when the Morse index of Mk{v) 
changes. 

For the polygonal equilibrium, we have decomposed the blocks Mq{i') and 
Mi{i') into moi(j^), mon(i^) and toii(j/), mi„(i/) respectively. Actually, 
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the orthogonal degree for the polygonal equilibrium has one element for each 
matrix mjk{i')- Thus, we may prove a more complete result because the group 
of symmetries is bigger. 

Definition 26 For j = 0, 1 and k — 1, ...^n, set nj}~(v) to he the Morse number 
of mjk{v), and define 

rj-jkiv) = o-{njk{v - p)~ njk{v + p)). 

The fixed point space of the isotropy group Fg x corresponds to the block 
'TiOn(O) = Bn- Since the generator of the kernel is Aid — T'„(— 77.^/^62), then, in 
this case, 62 is in the kernel of TOon(O). 

Let a be the sign of mon(O) in the space orthogonal to 62. Since /3„ = — a+si 
and an — a^si, then Bn ~ diag{{a + + si),0). Thus, for the polygonal 
equilibrium, a is the sign of {a + 1)(^ + si), |19^, p. 3220, 

a = (Tnip) = 1. 

In this case, we still need the hyperbolic condition near a bifurcation point. 
That is, we need to be sure that the matrices mofc(O) = B^ are invertible for 
— 1. In the paper [12], Proposition 23, we did prove that the block 
Bk is invertible except for one point fik € (si, 00). In fact, there is a bifurcation 
branch of relative equilibria from each pk for fc = l,...,n — 1. 

From the results of [2j, Remark 3.5, p. 259, we may state the following 
theorem for p ^ fii, pn-i- 

Theorem 27 If rjjkivk) is different from zero, then the polygonal equilibrium 
has a global bifurcation of periodic solutions from 2tt /vk, with isotropy group 
Zn{k) X Z2, for J = 0, and Zn{k) x Z2, for j = 1. 

7 Spectral analysis 

In the previous section, we have proved that there is a bifurcation of periodic 
solutions when the blocks of Mk{v) change their Morse index. The spatial 
and planar blocks were given in proposition (|10p . We will analyze the spatial 
spectrum for a general relative equilibrium, corresponding to Miiy). 

For the polygonal equilibrium, the blocks are given in propositions (|12p and 
(|16p . and the description of the isotropy groups was done in section four. In 
this section we analyze completely the spectrum of all these blocks mjk{i^)- 

7.1 Spatial spectrum for a general equilibrium 

The spatial block of a general relative equilibrium is 

were the matrix (a^ ) was given in proposition (jH)) 

The following result is well known in linear algebra. 
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Lemma 28 Let B{zQ,r) be the ball of center zq and radius r. Then, the 
spectrum a{A) of a matrix A = (aij)^_i is in the union ISi^i B{aii,ri) with 

The matrix {aij) has real eigenvalues because it is selfadjoint. And from 
the previous lemma, we have that the eigenvalues of the matrix (ay) are in the 
union [J"^^^ B{aii,\aii\), because an = — ^j^i0.ij. Since an is negative, then 
Mi(0) = {aij) must have n — 1 negative eigenvalues. 

Since the eigenvalues of Mi{v) are continuous, and Mi(0) has n — 1 negative 
eigenvalues, then the Morse number of Mi{v) satisfies ni{v) > n — 1 for any 
small p. Now, since ni(oo) = 0, then the matrix Mi(i/) must change its Morse 
index at n — 1 values of u. However, these n — 1 values could be the same, and 
we may assure the existence of only one point where the Morse index changes. 

Theorem 29 For the n-body problem, each relative equilibrium with a has 
at least one global bifurcation branch of periodic solutions with symmetries Z2. 
Generically, the equilibrium has n — I bifurcations of periodic eight solutions. 

In the particular case where all the bodies have the same mass, mj — m, 
let P be a matrix such that rn~^(ay ) = P~^AP, where A = diag{—v'l, —v^) 
with 

-vl < ... < <Pn = 0- 

Since 

Mi{iy) =mP-\v'^I + A)P, 

then the matrix Mi(z^) has eigenvalues \k = m{v'^ — u^). 

If all the life's arc different, then ni{vk ^ p) = k and ni{vk + p) = k — I. 
Thus, for k = 1, n — 1, the change of Morse index of Mi{i/) at the points 
is 

If the eigenvalue has multiplicity j, it is easy to prove that the change of 
the Morse index is r]i{i'k) = jo'. In any case, all bifurcation points have indices 
of the same sign, then for the equal mass n-body problem, these bifurcations 
branches must be non-admissible or connect to other relative equilibria. 

7.2 Spatial spectrum for the polygonal equilibrium 

Theorem 30 Define Vk = + for k S {1, ...,n — 1}, and Vn = ^ \i + n. 
Then, the polygonal equilibrium has a global bifurcation of periodic solutions 
from 2'jT/i'k with symmetries Zn{k) x Z2 for each k G {1, ...,n}. 

Proof. For k € {1, ...,n— 1}, we have that mikiiy) = — (t^ + Sk) changes only 
at Uk- Since u = 1, then r]i{i/k) = 1. For k = n, we have that 

det mi„(i/) = u^/j, [iP' — n)) 
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changes only at i^n- Since detTOi„ is negative for v < Vn, then nin{vn — p) = 1- 
Moreover, since ni„((X)) = 0, then ?7i„(i/„) — 1. The existence of the bifurcations 
follows from the bifurcation theorem. ■ 

Since the spatial bifurcations for the polygonal equilibrium have index rj = 1, 
then all the bifurcation branches must be non-admissible or go to the other 
equilibria. 

7.3 Planar spectrum for the polygonal equilibrium 

In the polygonal equilibrium a, the frequency is a; = fi + si. The equations have 
physical meaning for positive mass, /i > 0, although, the problem is well defined 
for fi > — si. In this section we will analyze the spectrum of the blocks mokiv), 
for fi > —si. 

Normalize the period of TOofe('^) by ^/uj and define mkiv) = mokiy^i^), then 

mfc(iy) = uj[i^'^I - 2iy{iJ)] + Bk for k G {2, n-2,n} and 

— uj[v'^diag{p.^ I) — 2i'diag{fi, ij)] + Bi. 

Notice that mkiy) is selfadjoint, so with real eigenvalues. Since mn-k{i^) = 
'fn-ki—i^), then the matrices TO„_fc(i^) and mk{—v) have the same spectrum. 
Thus, the Morse numbers nk{v) satisfy 

nn-k{v) = nk{-v). 

We cannot calculate explicitly the sums s^, but we shall use that the sums 
Sk are positive, satisfy Sk — Sn-k = Sn+k, and that Sk are increasing in k for 
k S {0, ...,n/2}. This last fact is proved in the appendix. The matrices Bk are 
given in the proposition ()12p . However, in order to simplify the computations, 
we shall restrict the analysis to the Newton case, that is with a = 2. 

7.3.1 Block k = n 

The block S„ is S„ (3/2)(/i + si){I + R) with R = diag{l, -1). Therefore, 

cr = sgn{e^Bnei) = 1, 

for /i > — si. 

Proposition 31 The matrix to„(j/) changes its Morse index only at the positive 
value u — 1 with 

Proof. The block m„(;/) is 

mn{v) = — 2u[iJ) + diag{3, 0)]. 

Thus, the determinant fi„(j/) = w^j/^ — 1) (i/ + 1) is zero only at ±1. Since 
dn{s) < 0, then n„(£) — 1. Moreover, since n„(oo) = 0, then 77„(1) = 1. ■ 

Therefore, there is a global branch of periodic solutions bifurcating from 
the polygonal equilibrium starting with the period 2tt and with symmetries 
Z„(n) X Za. 
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Remark 32 In fixed coordinates, the central body satisfies qoit) = and the 
other n bodies satisfy Uj+i{t) = e^^''Ui{t). Now, if v remains 1 on the branch, 
that is, if the frequency, of the solutions in fixed coordinates, isLO, then qj+i{t) 
• Thus, this bifurcation branch may be made of solutions moving on 

curves looking like ellipses. 

7.3.2 Blocks fee {2,..., n- 2} 

Define dk{v) to be the determinant of nik^i^). 

Proposition 33 The determinant of mkif) is 

= w^z^"^ + (2afe - w - si)u}v'^ - 4w7fei/ + + iibk, 

where 

ak = {si + - Pl-^l and bk = 3(si +ak+ fik) > 0. 

Proof. The block mk{v) is 

mk{y) = {ujv'^ + si + ak)I + (3m/2)(/ + i?) - /3fci? - {2vuj + 'yk){iJ)- 

Therefore, the determinant is 

dfe(z/) = (wz/^ + si+ak- Pk + 3/x)(a;zv^ + si + afe + Pk) 
- {2uu + 'ykf. 

Since to = ^ + si and dk{0) = bkfJ, + Uk, then 

dk{v) = uj'^v'^ + {2ak - oj - si)wi/^ - ioj^ki^ + ak + jJibk- 

■ 

We write the determinant in terms of uj as 

dk iy) = + 6u; — c, 

where 

a = u'^{i>^ - 1), 6 = !/^(2afe - si) - 4!/7/j + bk and c = sibk - ak- 

As a consequence of the next lemma, we shall prove that b^ + Aac is positive. 
Then v) is zero at exactly the two real solutions 



, , -b±VW+4ac 
"^("^ = 2-a • 

Lemma 34 For any v gR, we have that c and b{i') are positive and satisfy 

b^{v) > 9c. 
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Proof. Prom the definitions of a/j and bk we have 



Asibk = 3siSfe_i + 3siSfe+i - 6si + ISsiSk, 



and 



Attk = 2siSfe_i + 2siSk+i - 5sl - 9sl + ISsiSfe + Sk-iSk+i- 



Therefore, 



4c = 4(si6fc - ttk) = 9sfc - (sfe_i 



Sl) (Sfe+l - si) . 



Using the inequahty of the appendix 



(sfc+i - Sl) < 2sk - {sk-i 



Sl) < 2sk 



we conclude that 



7sl < 4c < 9s^ 



Since b^u) = 2i/(2afe — si)— 47^, then b'{i') = at i^o = 27fe/(2a;; — si). Since 
2ak — Sl > 0, then vq is a minimum of b{i'). Moreover, since {2ak — si) — 27/j = 
Sfe_i — Sl > 0, then z/q € (0, 1). Thus, 



From these inequalities one has that b'^{p) > 9(9s^/4) > 9c. ■ 

For i/^ > 1, since a is positive, then b"^ + Aac is positive. For € (0, 1), since 

a G (-1/4, 0) and b'^{v) > 4c, then 6^ + 4ac > 4c(l + a) > 0. Consequently, the 

two solutions co±{i') are real for 7^ —1, 0, 1. 

Remember that the equations of the bodies have meaning only for w = ^ + si 

> 0. For v'^ > 1, since 4ac is positive, then w+(i^) is positive and w-(z^) is 

negative. For i/"^ e (0, 1), since a and 4ac are negative, then the two roots 

uj±{i') have the sign of b. Since b is positive, then the two solutions w±(z^) are 

positive for i/^ e (0, 1). 

Defining /U±(i^) = oj±{i') — si, the matrix rukin,!^) changes its Morse index 

only at the two curves 



These curves satisfy the inequalities — ,si < /i+(i^) < M-(^^) for 1*^1 '= (O^l)- 
Moreover, since a — )■ when |:/| — )• {0, 1}, then /x_(i^) — )■ 00 when Iv] — > {0, 1}. 
Furthermore, since u)-\.{i') — )• when — )• 00, then /U+(i') — si when ^• 
00. 

Definition 35 Lei mo &e f/ie maximum of ^^{v) in R, m+ &e f/ie minimum of 
/x_(i/) (0,1), and m,- be the minimum of ^-{v) on (—1,0). 



b{vo) = M^lk - 47fc) + bk>bk- 27fc. 



Therefore, 



46(j/o) > 46ft - 87fc = 18sfc - 6si + 5sfc_i + s^+i > ISs^. 



/x+(z/) for 1/ G R, and /i_(i/) for i/ G (-1,0) U (0,1). 
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Figure 4: Graph of dfc(/i, u) = 0. 

Proposition 36 For k G {2,...,n— 2}, if n & (— si,/z/c), i/ie« mkif) changes 
its Morse index at a positive value Vk, with 

rjkivk) = 1- 

For k e {2, [n/2]}, if n {nk,rno), then wikiv) changes its Morse index 
at two positive values i>k <fk, with 

mi^k) = -1 and r^kifk) = 1- 

For k e {2, n — 2}, if fJ. G oo), then mk{v) changes its Morse index 
at two positive values Uk < t^k < ^, with 

r]k{vk) = -1 and r]k{i>k) = 1- 

Proof. Since /Xfc = —hk/ak is the only zero of dkifJ-jO), then ^k < mo- Notice 
that the function dfc(z^) + 4a;7fei^ is even in v, then 

1^) = dkifJ-, -v) - Su^ki^. 

Since -fk > for A: e [2,...,n/2) n N, then dk{n.,v) < dk{n.,-v) for v e R+. 
Thus, the maximum toq is reached for some > 0, and m+ < to_. 

Since mk{i^) is a 2 x 2 matrix, then the Morse number is nk{v) = 1 when 
rffc(z^) is negative. Hence, the Morse index is nk{n, = 1 in f2, where 

Since nfc(i/) = for large enough, and since mk{v) changes its Morse index 
only at dVl, then j/) = in f2^. 

As we have seen before, cr = 1 for /x > — si, then, for the first case, we have 
that rjk{vk) = 1 — 0. For the other two cases, there are two values, Vk < I'k, 
with r?fe(i^fe) =0-1 and %(z^fe) = 1-0. ■ 
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Remark 37 Since 

d'l{i>) = 2{6J^v'^ + {2ak - H- 2si)uj), 

then d'l^{v) is positive for ii, < 2{ak — si). Thus, for /i G {—si,2{ak — si)) the 
matrix mk{v) changes its Morse index only at u^, for the first case. Now, since 
dkif) is a polynomial in v of degree 4, then dk{y) has at most four zeros. This 
is the case for ^ > max{m+, m_}, and then the two positive solutions for the 
third case are the only ones. 

From the bifurcation theorem we have the foUowing result, for n > 4. 

Theorem 38 The polygonal equilibrium has a global bifurcation of planar pe- 
riodic solutions with symmetries Zn{k) for each k e {2,...,n — 2} and fi G 
(— sij^fe)- When /i G {—si,2{a.k — si)), these bifurcation branches are non- 
admissible or go to another equilibrium . 

For each k G {2, [n/2]} and /x e (/Ufc, mo), and for each k G {2, ...,n— 2} 
and II G (m^, oo), the polygonal equilibrium has two global bifurcations of planar 
periodic solutions with symmetries Z„(fc). 

Remark 39 // we had proved that there are only three solutions for dk {v) = 
and d'^ {v) = 0, then these solutions would have been mo and m± . In that case, 
the only points where the Morse index changes are those of the previous theorem 

We can prove this result for k = n/2, and large n, since 7„/2 = 0. In fact, 
one has that m_|_ = m_, and mo = fJ'n/2 '■s large enough. Since 7„/2 = 0, the 
graph for dn/2 is symmetric with respect to the fi-axis and one has two curves 
for n as a function of v. 

In order to compute m+ = m,^ , notice that 

d!^{v) = 2i/a;2(2!/^ + w"^(2a„/2 - 2si - /x)) = 0. 

Substituting = — a;~^(2a„/2 — 2si — /u)/2 in dkiy) = 0, we obtain 

(2a„/2 - 2si - iif - 4(a„/2 + /U^„/2) = 0. 

From this quadratic equation, we get that m+ = 6 + \fW^Vc, with 

b = 2(6„/2 + an/2 - si) and c = 4(a„/2 - (a„/2 - si)^). 

The other root, m = b — \Jb + c will be valid provided the expression for 
is positive, giving at most two values for ji, and v > 0, where the derivative of 
ji with respect to v is fi' = 0. At these points one has that 

IJ," = -dvv/dui = 8w(2a„/2 - 2si - li)/{b - jj). 

In particular, n" > at m+ and ^" < at rh. 

It is easy to see that rh < 2(q!„/2 ^ si) (and there the corresponding is 
negative) if and only if iJin/2 satisfies the same inequality and, in this case, iJin/2 
is the only maximun of this piece of the graph. 

For n large, the asymptotics of the appendix show that /i„/2 behaves like 
or? 1 2), while 2(a„/2~si) *s ^^^e an'^ , which implies that, forn large, mo = A*„/2- 
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7.3.3 Blocks fc e {1, n - 1} 

From the definition, the block toi(j^) is 

/ n{u:v'^ ~2vuj + si+ n + nl2) -2{n/2f''^ ^ -{n/2f''^ni 

—2{n/2)^^'^^ oji'^ + si + ai + 3/i aii + 2vuji 

y /ii ^aii — 2vuji uji>'^ + si + ai 

Let us define ai and bi as 

ai = (2si + n)(2si + S2)/4 and &i = 3(4si + S2 - 2n)/4. 

The determinant of rni(z^) may be written as 

det mi(zy) — LUji^v — 

where di{iy) is the polynomial 

= - 1)^2 + [(S2 - 2si + n) iy'^/2 - (sa - n)iy] a; + (oi + 

Therefore, we need to find the points where di{h') is zero. 

Lemma 40 The functions bi, 2ai — si, si — n, S2 — n and 8(si — n) + 9s2 are 
positive, at least for n > 1071. For n > 3, we get numerically that si — n is 
negative if n < 472 and positive if n > 473, that S2 — n is negative if n < 11 
and positive if n > 12, and 



sgn{bi) — sgn{2ai — si) — sgn{8{si — n) + 9s2) 



— 1 for ri < 6 
+1 for n >7 



Proof. Since sin x < x, then 

1 1 1 1 1 , . ln(n/2) 

si = -) ^>- > > - n n/2) ^ ' ' . 

4 ^ sin jC/2 ~ 2 ^ jC 2 - ^ ' ' 2tt 

Since ln(n/2) > 2tt for n > 1071, then si > n at least for n > 1071. 

Using that §2 > si and bi = 4(2(si — 7i) + S2 + 2si)/4, we conclude that bi, 
82 — n and 8(si — n) + 9s2 are positive for n > 1071. Finally, from proposition 
([57]). we have that S2 = 4si — si with si = sin(jC/2) < fi, then, 

2ai — si — S2/2 — si = si — si/2 > si — n > 

for n > 1071. ■ 

We shall analyze the graph of di{i') = only for positive masses, fi > 0. 

Proposition 41 For n > 7, the polynomial di {11, v) is zero only for the function 

/ioH : (-1,0)U(0,1)^K+, 

where iJ-oii') ^ 00 when 0, 1 . 

For n e {3, 4, 5, 6}, the polynomial i') is zero only on the curve 

(Aici'o) X (-c^,l), 

where /io(i) — >■ 00 and i^oit) — > ±1 when t — )• ±00. 
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(a) n > 6 (b) n = 3, 4, 5, 6 

Figure 5: Graph of v) = 0. 



Proof. Since uj = fi + Si, 

di{n, v) = {iP- - -V {s\h\ - a\ -V c\ -V s\ [y^ - !y^)]/x/si + ci, 
where Si6i — ai = (8si + S2) (si — n) /4 and 

ci(i.) = [2si(i.- l)2+n][2si(z. + l)2 + .s2]/4. 

For each fixed i^, set di{i') be the value of the determinant. Then, its deriva- 
tive d[ {fj,) is zero at the only critical point fic such that 

— (si6i -ai+ci+sl (1^^ - iy^))fj,c = -2i/^ - 1) i4- 



For < 1, since didic,!^) 



1) fi"^ + ci is positive, then the polyno- 



mial has at most one zero for fj, > fj,c- Since di/u 



v'^{v'^ — 1) when 



H ^ QO, and di{Q,v) = ci is positive, then the polynomial di has exactly one 
zero at /xo(j^) > for each G (0, 1). 

In order to obtain the limits of the function Hoiv), when — )■ 0, 1, we need 
to study the zeros of di{ix,v) for \v\ = 0, 1. For v = the polynomial 0) 
is zero at /Uq = — ai/61. Since ai is positive and 61 is negative for n < 6 only, 
then yUo is positive for n < 6 only. For 1^ = 1, the polynomial di{ii, 1) is zero at 
the point /U+i = —n. For v = —1, the polynomial di{fi, —1) is zero at 

M-i = -S2(8si + n)l (8si - 8n -h 9s2), 

which is positive for n < 6 only. 

Hence, for n>l, the polynomial di(pL,v) is positive at \v\ = 0,1, for any 
/X > 0. Therefore, the functions Hoiy) — > 00 when — )■ 0,1. Moreover, since 
ci{v) is increasing for v > 1, then di(i^) > di(l) > 0, for v > 1. Similarly, one 
has that di{v) > 0, for < — 1. We conclude that the only zeros of v) for 
/X > are for the function Ho^iy). 
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For n < 6, we may prove, in a similar way, that di{fi, i>) is positive, for v > 1. 
As a consequence, the function /Uo(i^) — >■ oo when i/ — >■ 1. Also, the function 
/io(i') has a continuous extension, to u = 0,-1, by /io(0) = A*o and no{—'^) = 
Furthermore, the continuum of zeros of j^) at 1) crosses the 

line (/X, — 1) only once, then this continuum must go to (/U, i^) = (oo,— 1). This 
statement follows from the fact that di{fi, u) is positive for i/ < — 1 and = 0, 
for /X > and \v\ large, and iov v < —1 and /x large. ■ 

Remark 42 T/ie equation di{jjL,v) = is quadratic in fx, thus we can find 
explicitly two roots /i±(z/)o/ i^) = 0. In principle, we could analyze these 
roots analogously to the case k G {2,...,n— 2}. However, there are many cases 
depending on n. The main case is for n > 473, where we may prove that the 
graph of d\{fi, v) =Q looks like the graph of dk{i^, v) = 0, for ji > —s\. 

Definition 43 For n > 7, define m+ to be the minimum, of the function ^H]{v) 
in (0,1), and m_ the minimum in (—1,0). For n € {3,..., 6} define mo to be 
the minimum of /io(f) for t e M. 

Proposition 44 For masses /i > we have the following cases: 

(a) For n € {3,..., 6}, the matrix mi{u) changes its Morse index at the two 

values v- < D- < if iJ. £ {mo,i^i), and at u- < Q < uj^ < 1 if ji £ 
(Mi,oo). 

(b) For n > 12, the matrix mi(z/) changes its Morse index at < 1/+ < 1/+ < 1 

if H & (m+,m_), and at the four values of v 

-1< v- < D- <0 < 9+ < 1^+ < 1 

if IJ: & (m_,oo). For n G {7, ...,11}, this statement is true with m+ and 
rri- interchanged. 

In both cases, we have that 

mi'^±) = -1 and ?7i(j/±) = 1. 

Proof. For n > 12, using S2 > n, we may prove that ci(i^) < ci(— i/). Then, 
rfi(/Lt, v) < rfi(/i, —v), for v > 0. Therefore, m+ < m_, for n > 12. In a similar 
way. wc may prove that m_ > m+, for 7 < n < 11. 
The trace of mi (i^) is 

Ti = fiujiu - If + 2w(z/2 + 1) + /x(n/2 + 1) + 2ai. 

Since Ti is positive for /i > 0, then the Morse index satisfies ni(r^) = 1, when 
detmi(!^) is negative, and ni & {0,2}, when detmi(i/) is positive. 

Let fl be the connected component where rfi(/x, v) is negative, for n > 7 this 
set is = {(iiji/) : (x > iio{i^)}- From the definition, we have that ui = 1 in 
ri — {(/i, ±1)}. Now, for /X = 0, the matrix mi(0, ly) has eigenvalues 0, si {ly — 1)^ 
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and si (j^ + 1)^ + 2ai, hence ni{e, v) < 1, for small e. Since ni € {0, 2} and 
ni < 1 in — ±1)}, then m = there. Thus, we have that 771 (±1) = — 0, 
r]i{i'±) = 0—1 and ?7i(i'±) = 1 — 0. ■ 

Prom the bifurcation theorem, we get the following result: 

Theorem 45 For n> 12, if ^> m+, the polygonal equilibrium has two global 

bifurcation branches of planar periodic solutions with symmetries Z„(l). If fJ. > 
rri-, there are two more bifurcation branches with symmetries Z„(n — 1). For 
n G {7, .... 11}, the same statement is true interchanging to_|_ and m_ . 

For n G {3, 6}, if 11 > /ii, the polygonal equilibrium has one global bifurca- 
tion branch of planar periodic solutions with symmetries Z„(l), and one more 
with symmetries Z„(n — 1). If H & (mo,/xi), there are two global bifurcation 
branches of planar periodic solutions with symmetries Zn{n— 1). 

Remark 46 When /U = 0, the central body does not exist, then the block mi {u) 
is 

I \ ( sxv^ -\- s\ -\- 0L\ ci\i^2vs\i \ 
"mxiv) = . _ . 2 , 

^ ' y —Q.\i — 2vs\i Sif + Si + ai J 

In this case mi(i^) has eigenvalues si{i' — 1)^ and Si(z^ + 1)^ + 2ai. Then 
ni {v) = for all v, and there is no bifurcation in this case. 



8 Remarks 

8.1 The case n = 2 

The irreducible representations for n = 2 are different from those for larger n, 
due to the action of Z2. This fact gives a change in the planar representation 
for k = 1, but the change of variables for fc = 2 remains the same as for larger 
n's. In particular, the bifurcation results for the spatial case and for the planar 
case, with k = 2 are similar to the ones already given. 
For k = 1, define the isomorphism Ti : — )• Wi as 

Ti{v,w) = (w,2-^/^u;,2-^/^u;) with 
Wi = {{v, w,w) : v,w G 

It is not difficult to find the action of Z2 on Wi and to compute the Hessian 
of the potential V. The matrix B2 is the same as before, but Bi is now a 4 x 4 
complex matrix. For a = 2, one has that Si = 1/4 and u) = n + l/A. The matrix 
mi{i^) is 



/X (uii^'^ +11 + 17/4) 
2ih'ujfj, 
-2(n/2)'/V 




IjOV 



—2iuu)fi 


(2)1/2 



7/4) 



-2(2)1/^ 


' +3/i4 
2ivu! 



1/4 



\ 

(2) 1/2/x 

—2ivLU 

2^1/4 J 
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The determinant of mi(z/) may be written as 

detmi(iy) = fi^ {Afi + if {v^ - lfdi{i^), 
where di(i^) is the polynomial 

di{iy) = {Wfi^ + 8fi+l)iy^+ (-IG/i^ + 20fi + 6) - (84/i + 119) . 
For fj, > 0, di^iy) is zero only for ivi, where i^i is the positive solution of 

lyf = (2/1 - 3 + 2 (/i^ + 18m + 32)^''^) /(4/i + 1). 

Since > 1 for /i > 0, it is enough to know the eigenvalues of toi(j^) 
for one value of fi, for instance, jj, = 1. These eigenvalues are 5{h' ± 1)^ and 
(5i/2± 2(251/2 + 81)1/2 + 11) /4. 

Thus, one sees that ni = 1 for —i^i < v < vi and ni = otherwise. Hence, 
Vii^i) ~ 1 and 771(1) = 0. Therefore, 

Theorem 47 There is a global bifurcation branch of planar periodic solutions, 
from 2t:/vi, with symmetries ^2(1). There is also a global bifurcation of planar 
periodic solutions, from the period 2tt, and symmetry ^2(2) x Z2; that is 

UQ{t) = and U2(t) — —Ui{t). 



Finally, there are two spatial bifurcation branches, from v = -^//i + 1/4 and 
V = y'/x + 2, with symmetries ^2(1) x Z2 and ^2(2) x Z2, respectively, that is, 
the first type is such that uo(t) = and zoit) = 0, 

U2{t) — —ui{t) and Z2{t) — —ziit), 

while, the second type is with uo{t) = and zo{t + vr) = zo{t), 

U2(t) = —ui{t) — —ui{t + tt) and Z2{t) = zi{t) — —zi{t + tt). 



8.2 Resonances 

We have proved that there is a bifurcation of periodic solutions with the sym- 
metries (O when the block Mi{i/) changes its Morse index at vi. However, we 
cannot deduce, from the symmetries ([7]), that the solutions are truly spatial, 
because all the spatial coordinates can be zero: Zj{t) = 0. In fact, there is a 
non-resonance condition to assure that the solutions are truly spatial. 

Proposition 48 If AIi{i>) changes its Morse index at vi, and the matrix Mq{21i'i) 
is invertible for all I S N, then the equilibrium xq has a bifurcation of truly spa- 
tial periodic solutions from vi with symmetries ^ . 
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Proof. Wc need to prove that near (a, Vk) there are no planar solutions with 
symmetry Z2, i.e. 

Xj{t) = Xj{t + tt), yj{t) = yj{t + tt) and Zj{t) = 0. 

These solutions have isotropy group Z2 x Z2, one Z2 is generated by k and the 
other by TT S 5^. 

Let f'{a) be the derivative of the bifurcation operator f{x), at a, and define 
/'(a)|yyH as the derivative in the fixed point space of = Z2 x Z2. The linear 
map /'(a)|w" has a block of the form Mq{21v), for each Fourier mode I. 

By assumption, the blocks Mo(2/i/i) are invertible, then the derivative /'(a)|w« 
is invertible. Using the implicit function theorem, we may conclude the non- 
existence of periodic solutions with isotropy group Z2 x Z2 near (a^vi). As a 
consequence, the bifurcation from vi must be spatial. ■ 

For the polygonal equilibrium, we have proved before that the block mik{vk) 
changes its Morse index at I'k = \J + Sfe- Restricting /'(a) to the group Z„(fc) x 
Z2 x Z2, we may prove, analogously to the general case, that the bifurcation 
from Vk is truly spatial if mnk^lvk) is invertible for all Z's. 

For k G {1, ...,n — 1}, we have proved that the block mofe(z^) is invertible 
for \v\ > ^[i + si, when n > mo: = 1 in the previous graph corresponds 
\,o V = \J s\ before the normalization used in the study of mk{v). Since 
Hvk > \/ \i + si for all I's, then mok{2li/k) is invertible for all I's, when /x > mo. 
Therefore, the bifurcation from Vk is truly spatial for k G {l,...,n — 1}, when 
^ > mo. 

If one uses the same normalization done for the planar spectrum, that is 
'mokiV^'^)^ for the spatial spectrum, then one has that, for k = 1, n — 1, 

mik{\/ujv) = i^'^ifJ, + si) - (m + Sk), 

and the corresponding expression for /c = n. 

Then, the block mik changes its Morse index at Uk = + Sfc)/(/U + si). 
Restricting /'(a) to the group Z„(fc) x Z2 x Z2, we may prove, analogously to 
the general case, that the bifurcation from Uk is truly spatial if mofc(2/i//c) is 
invertible for all Vs. 

For k G {1, .... n — 1}, we have proved that the block mofe(j^) is invertible for 

> 1, when n > mo. Since > 1 for all /'s, then mofe(2lz/fe) is invertible 
for all Vs, when n > mo. Therefore, the bifurcation from Uk is truly spatial for 
k G {l,...,n — 1}, when /i > mo. For other values of /i < mo, one has that 
^fc = (a* + + si), and w = /x + si = {sk — si)/(vk'^ — 1). The equation 

dki^lvk) = leads to finding the zeros of a 4th degree polynomial in Vk, with 
coefficients which depend on I. For the piece of the graph where I'o < Ij one 
cannot have the resonance vo — "^-Ivk and, for the rest of the graph which is 
asymptotic to the line n = — Si, there will be only a finite number of Z's, where 
one could have the equality of the frequencies. As a consequence, there is only 
a finite number of possible resonances, for = — 1. 

For k = n, v^^ = (y^ + n)l{\x + si) and vq = 1, so that there is at most a 
finite number of Z's and positive /x's which satisfy the relation vo = Hvn- 
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Remark 49 // one has a resonance between spatial frequencies, for ki and k2, 
such that V2 — Ivi, then one has the relation — 1) = — I'^Ski, where 3^2 
is replaced by n for k2 = n. If I = 1 and k2 ^ n, then k^ — n — ki, and the 
bifurcating orbit for k2 is the one for ki taken in the opposite direction, due to 
the symmetry k. 

If I — 1 and /c2 = n, then one needs Ski — n for a 1 — 1 resonance. We 
know that 3^/2 > ••• > ^i, and we have seen, in the last lemma , that si > n 
for n > 472, that 32 > n > si for 12 < n < 472 and n > 32 for 3 < 71 < 11. 
Furthermore, numerically one computes n > S3 for 6 < n < 7 and S3 > n for 
8 < n < 11. Since all these inequalities between sq = n and Sk are strict, then 
we may conclude that there are strict inequalities between vq OLi^d Vki , that is 
there is no 1 — 1 resonance in this case. 

For k2 < n, we may assume that I < ki < k2 < n/2 and we need, for /i > 0, 
to have P < s^^l^ki- Using the recurrence formulae of the appendix, we see 
that Sfej/sfcj^ < fc| — fc^ + 1 < fej. Hence, the number of possible subharmonic 
resonances is finite, as well as the number of possible solutions to the relation 

tl={sk2-P3k,)/{P-l). 

Remark 50 For resonances between planar frequencies, if V2 = Ivi anddkj{i'j) = 
for j = 1, 2, it is easy to see that v on the graph ofdkii') — is bounded above, 
for i-J, > and away from for all positive bounded fi 's, except for small neigh- 
borhoods of Ilk- Thus, the number of possible I's is finite. Now, if kj < n, 
dk2 (f^i) ~ l^dki (^1) polynomial of degree 2 m /i and vi, which can be solved 
for /i in terms of vi and dk^{vi) gives a polynomial of degree 16 in vi, if I > I, 
and of degree A if I = 1. In any case, we obtain at most 16 solutions for vi 
and 32 for ^i, i.e., a finite number of possible resonances, with I bounded if jj, is 
positive, bounded, and not in a neighborhood of , where vi is small. 

Note that, if k2 — n, then 1^2 = 1 and there is a resonance for vi — l/l and 
the corresponding /i, with I going to infinity as /i goes to fiki- For ki < n/2, 
assume that the graph for dj.^ (v) covers a small interval [/ij.-^ , fi] . Then, if I is 
such that Ivki (A) ^ riUk2 (A), the graphs for dki {Iv) will cut the graph for dk2 {v), 
for any I > I, in the interval, at fi{l) with the limit ^ki, giving rise to an infinite 
number of resonances. This situation will happen if fik2 > Mfci , which is the case 
for n large enough (see remark 29 in fl9f). But, for fi different from iiki, one 
has a finite number of possible resonances. 

A similar argument may be done for ki > n/2 and intervals at the left of 

Note that, at fiki one has a two-dimensional kernel for the linearization at 
the relative equilibrium (for ki and n ) and the orbit is no longer hyperbolic, so 
one needs a different computation of the bifurcation index. 

Remark 51 Assume that there is no resonance, that is that dj{lvk) 7^ 0, for 
j k and any I > 1, and Vk ^ 1. Then the equations for j ^ k, j ^ n 
may be solved close to the bifurcation orbit, by an implicit function argument in 
terms of the component for k. The equation for n, which has a one dimensional 
kernel, is solved by a Poincare section argument which is used in f24f in order 
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to compute the orthogonal index for an isolated hyperbolic orbit (theorem 3.1 p. 
247 in the above reference). One has then a reduction to one mode, a single k 
and a one- dimensional complex orthogonal equation. The non-trivial solutions 
of this reduced system have a least period equal to 27r, which implies that, except 
for a finite number of /i 's, if different from fik , the branches which bifurcate 
from the relative equilibrium, for fixed /i, will be different. 

8.3 Linear stability 

The linearization of the equation at the equiUbrium xq, is 

Mil + 2y/i:JMJu = D^V{xa)u. (12) 

For this problem, an equilibrium is linearly stable if the linear equation p2)) has 
6(n + 1) linearly independent solutions of the form u{t) = e^^^'wo with i'q S M, 
for n + 1 bodies. 

Now, if I'o G is a bifurcation value, that is the determinant is at z^qi then 
there is a vector vq such that Af(i'o)^o = 0. Since 

-Mil - + D'^V{xq)u = e'''°*M{i^o)vo = 0, 

then the function u{t) = e^'^°*vo is a solution of the linearized equation p2|) . 

Therefore, for each point vq where Miy) is singular, the equilibrium has a 
direction where the linearized equation is stable. Now, vq may be a multiple 
zero of the characteristic determinant and one could have a polynomial increase 
in t. A weaker definition of stability is that of spectral stability, where the 
equilibrium xq is spectrally stable if M{v) is non-invertible at Q{n + 1) values 
of V, counted with multiplicities. 

As a consequence, we may deduce the spectral stability from the spectral 
analysis which we did to prove bifurcation. We shall recover, with somewhat 
different arguments, the results of [31], of [3S] and others. 

For the polygonal equilibrium, we may give a full description of stability. 
We have proved that each spatial block mi]^{v) is non-invertible at exactly one 
positive value, for k S and mi„(0) has a one-dimensional kernel, 

with a double zero. Then, the spatial spectrum gives n stable solutions. The 
kernel of mi„(0) is generated by the vector {\,y/n), which gives, under the 
transformation P the vector (1, 1, , , , 1) = e. The eigenvectors corresponding to 
the double root are {a-\-bt)e. In fact, since = 0, the vector e generates the 

kernel of the vertical Hessian. Note that, if one sums all the vertical components 
of the equations, one gets '^z'j — Q and, if the initial conditions for ^ Zj are 
0, for position and speed, then this sum remains 0. Hence, by imposing these 
restrictions, one obtains a linear stability in the vertical direction. 

Next, we analyze the planar spectrum. 

Proposition 52 Let nit, be the largest of the m_|- 's, for k = 1, n. Then, the 
polygonal equilibrium is spectrally stable only if n > 7 and 11 > . It will be 
linearly stable if one restricts vertical translations, rotations around the center 
of mass and scaling of the frequency. 
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Proof. For k S {2, ...,n — 2}, we have proved that the 2x2 matrix mk{v) is 
non-invertible at four values of when /i > m* . For k = n, the 2x2 matrix 
rn„(i^) is non-invertible at one positive and one negative value, and is a double 
root, with a kernel generated by the vector (0, 1). Applying the transformation 
P, one gets the eigenvectors Jxq and 2xo — 3^/uJtJxo. Now, consider the relation 
\7V{uja~^ ,axo) = 0, for any positive a, and differentiate it, with respect to a, 
in order to obtain that D'^V{xo)xo = SMujxq. Since Jxq generates the kernel 
of D'^V{xo), one verifies that the two vectors are solutions of the linearized 
equation. Hence, the behavior at = is a consequence of the symmetry of 
rotations in the (x, ?/)-plane and scaling. Thus, this does not destabilize the 
polygonal equilibrium, if one asks for this natural restriction. 

For k E {l,n— 1}, we proved that the 3x3 matrix mfc(i/) is non-invcrtible 
at three values, with = 1 as a double root, for n G {3, 6}. Hence, since the 
total count of zeros is An, the polygonal equilibrium is never spectrally stable for 
n e {3, 6}. Similarly for n — 2, the 4x4 matrix toi(i^) has a characteristic 
determinant with six real and two purely imaginary zeros, giving an exponential 
instability. 

For n > 7, we proved that mi{i>) is non-invertible at six values of v, with 
V = 1 being a double root, ii fi > m^,. Therefore, the polygonal equilibrium is 
spectrally stable only if /x > m*, since one has the complete count of 4(n -I- 1) 
roots for the planar spectrum. 

Now, it is easy to prove that, for v ^ \, the matrix mi(l) has a kernel 
generated by the vector (-\/2/n, 1, i), which gives, under the transformation 
P, the vector v = l/A/n((l, z), (1, i), (1, i)). It is then easy to prove that 
(a -|- ht)e^*v are solutions of the linearized equations. In fact, if one takes the 
sum of all the 2x2 matrices X^^- Aij — 0, one has that v is in the kernel of the 
Hessian of the Newton's potential, and (a + 6t)e\/"*(l, i) generates the four (in 
its real and imaginary parts) solutions of the equation 

X + 2\fujJx — ux — 0. 

Furthermore, since Newton's potential is invariant under rotations in the plane 
of all the masses, then its gradient is equi variant and its Hessian has v as 
generator of its kernel. This explains the restriction on the rotations. ■ 

Assume the assertion, in the paper [32], that the largest of the 771+ 's corre- 
sponds to the block rriniiiy). In the remark ((39|). we have found, for fc = n/2, 
that m+ = 6 + + c with 

h = 2(fe„/2 + q;„/2 - si) and c = 4(a„/2 - (a„/2 - si)^). 

Next, we wish to find the estimate of to+, given in [32], on the stability of 
Saturn's rings. 

Lemma 53 [32]: the sums Sk satisfy the limits 

..Si 1 Sn/2 1 1 

iim — = — ana lim — — — a — — - > — -7. 

n^oo nlnrt 2tt n^oo n"^ 2n^ {2k — ly 
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From proposition ([57)) . given in the appendix, where si is defined and a is 
taken to be 2, we have 

n/2-l 

Sn/2 - Sn/2-1 ^ {n ~ l)si - ^ Si- 

1=1 

Using approximations by integrals, we can prove that si/n is finite and that 
si/n — oo when n — >■ oo. Since s; = 0(si), when n — )■ oo, we have that 

Sn/2 - Sn/2-i = 0{nsi) = 0(n^ Inn) = o(n^). 

Consequently, we have that lim„_i.oo = 
Using these limits and ^n/2 = 0, we have that 

an/2/n^ -> 0-/2 and Pn/2/n^ (3/2)cr, 

when n — > 00. From the definitions of a„/2 and 6„/2, we obtain that a„/2/n® — >■ 
— 2ct^ and fo„/2/f^'^ ^ 6cr. Therefore, fo/n'^ — > ISct and c/n^ — > — 9cr^, when 
n — 00. We conclude that 

lim m+/n^ = (13 + 4%/lO)cr. 

n— f 00 

This limit is the one found in [32] in order to estimate the stability of Saturn's 
rings. 



8.4 The problem of n-charges 

We wish to analyze the movement of n particles with charge —1 interacting 
with a fixed nucleus with charge g > 0. We suppose that the gravitational 
forces are much smaller than Coulomb's forces. Thus, these equations may be 
a classical model for the atom. Since electrons and protons have equal charge 
with different signs, then q = n for a non-ionized atom. 

Let Xj be the positions of the negative charges for j £ {l,...,n}, then the 
equations describing the movement of the charges are 

n 

.. ^ J > J Xi 

where the first term represents the interaction with the fixed nucleus. 
In rotating coordinates, Xj{t) — e^*'^Uj{t), the equations become 

ilj + 2^/ujJuj — VV^(m) with 

F(u) = ^ V ||/z.,ir - V j-^—r + E A- 
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Let u be {ui, ...,u„), then 

u + 2^^; = Vt/(u). 

We wish to show the similarities of this problem with the (n + l)-body 
problem. If we set wq = in the potential of the (n + 1) bodies, the potential 
for the bodies is 

We have proved before that the polygonal equilibrium (0, ai, a„), with aj = 
e'-''', is a critical point of the potential 1^(0,^). 

If we choose w = g — Si, in V{u), and put /x = — g in V{0,u), then the 
potential for the n-charges satisfies 

V{u) = -V{0,u). 

From this equality, we have that a = (ai, a„) is a critical point of V{u), with 
oj = q — Si > 0, because a is a critical point of V{0, u). 

Thus, we may perform a similar analysis to the bifurcation of periodic solu- 
tions from the polygonal equilibrium in the body problem. The only difference 
is in the spectrum. Hence, we shall focus on the spectrum of the polygonal 
equilibrium for the charges. 

The block corresponding to the Fourier modes of the charges is 

M(iy) = _ 2i'y/U{iJ) + D^V{a). 

Using the computation of D'^V and the equality D'^V{u;q) = —D'^V{0,u; —q), 
and after the change of variables, we obtain the decomposition of M{i^) into two 
blocks. Next, the planar spectrum is decomposed into the blocks 

mofe(i^) = i^^/ - 2i^y/u}{iJ) - Bki-q), 

for k e {1, ...,n}, where Bk{^) are the matrices given in the body problem 

Bk{n) = (3/2)(7 + R)^i + (si + ak)I - PkR - IkU- 

While the spatial spectrum is decomposed into the blocks 

■mik{v) = i/^ + (-g + Sfe) 

for k e {1, n}. 

For k e {1, ■■■,n}, \i q> Sk, then the block mife(z^) changes its Morse index 
at the value 

j^ife = Vq- Sk 

with Tjikivik) = 1- 
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Proposition 54 For each k G {l,...,n}, the block mok{i^) changes its Morse 

index at a positive value Vk{lA for q G (si, oo) with 

Proof. For k = n, wc have that detmon('^) = — 1), thus detmon(z^) is 

zero only at = 1. For k G {1, ...,n — 1}, the determinant of mofc(i^) is 

dkil, ^) = w^!/'' - (2afe + w - sijwv^ + 4w7fc!/ + - g6fc 

with u> = q — si. Since dk{q, 0) = — qhk is negative and dk{v) is positive for 
large i^, then there is a positive value Vk where mofe(At) changes its Morse index. 
For 5 > si, we have cr = 1 since mon(O) = diag{3{q — si), 0). 

Since dk{q,0) = ak — qbk is negative, then nfc(O) = 1, and since n(oo) = 0, 
then rjkii^k) = 1-0. ■ 

Theorem 55 For n>2 and q G (si, oo), the polygonal equilibrium has a global 
bifurcation of planar periodic solutions with isotropy group Z„ (fc) x Z2 for each 
k G {1, ...,n}. If q > Sk, there is a global bifurcation of spatial periodic solutions 
with isotropy group Z„(fc) x Z2 for each k G {1, n} and starting at vik- 

Remark 56 It is possible prove that the function Vkil^) decreases from, Vk{si) = 
00 to Vk{oo) — 1, and that Vk is unique for k G {2, ...,n — 2,n}, and for k G 
{l,n— 1} if n > 7. In those cases the bifurcation is non- admissible or goes to 
another equilibrium. Note also that, if q = n, then one may verify numerically 
that si < q only for n < 473. 

9 Appendix 

Let ^ = 27r/n, for a fixed a G [l,oo), we define the sums Sk and Sk as 

sin^(fcjC/2) , _ sin^(fcjC/2) 
- 2« A. sin"+i(iC/2) - 2-^ sin"-(iC/2) ' 

Then, the sums Sk are always positive and satisfy 

^k ~ ^n+k ~ ^n—k- 

Proposition 57 The sums Sk satisfy the following three recurrence formulae 

Sk+i - '^Sk + Sk-1 = 2si - Sk 
fc-i 

Sk = k^si — Isk-i and 
1=1 

k 

Sk+l -Sk = {2k + l)si - ^ . 

1=1 
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Proof. Write Sk as 



1 1 — cos{kjC,) 



^-sin"-i(jC/2) l-cos(jC) ■ 
Using geometric series, we obtain that 



1 - cos(jC) 1 - e^i'^ 1 - e-'i'^ ^ ^ 
Therefore, the difference of the sums satisfies 



e"- 

;=o m=0 



Since the sum of exponents is 

h, h h k 

^ e^^f^i^ ^ cosjX= 5] (l-2sin2(j/iC/2)) = (2fc+l)-4^sin2(jK/2), 

h— — k h — — k h— — k h—1 

then, 

h=l j=l h=l 

Iterating the previous formula we obtain the other two recurrence formulae. ■ 
We have used the results about the s^'s for the vortex problem, a = 1, in 
our paper [20j . Now, we may conclude the following corollary. 

Corollary 58 Using the fact that Sk and 4sfc — Sk are positive, one has 

Sk+i - 2sk + Sk-i < 2si and Sk+i + 2sk + Sfc-i > 2si. 

For the n-body problem, a = 2, we cannot find explicitly the sums Sk- 
However, we may prove that the sums Sk are increasing, for k £ {0, [n/2]}. 

Lemma 59 Let s be the sum s — — J2j=o shi(A — l/2)jCj then s — — cot(fc — 
l/2)C/2. 

Proof. Let ^ = {k - l/2)(, then 



= -^sin,e = -l^(e--e 



j=0 3=0 



n~l 

2i 
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Using geometric series, we have 

1 / (1 - e*«") (l - e-*«") \ 
^ ~ 'Yi \ {\- e*«) (1 - e-»«) j 



Thus, 

sin^ — sinn^ + sin(n — 1)^ 
* ~ 2(1- cos ■ 

Since — 27r(fc— 1/2) = tt modulus 27r, then sinn^ = and sin(n— 1)^ = sinf . 
Hence, 

sin^ 



1 — cos ^ 



■cot(C/2). 



Proposition 60 The sums Sk for k G {1, [»t./2]} are increasing, 

Sk > Sk-l- 

Proof. Define dk, d'j, and dj! as 

dk = Sk+i - Sk, d[ ^ dk - dk-i and d'l = d'^ - d[_j^. 

The goal is to prove that d'^ are negative, and that dk are positive, for k G 
{l,...,[n/2]}. 

From the equality (fT^ . we have, for o; = 2, that 

^ n— 1 k 

' 4 A. 3inO-C/2) • 
Therefore, the differences d'^. are 



, _ 1 2cosjfcC 



4^sin(jC/2)- 



Thus, 

,„ ^ 1 cos jkC - cos j{k - 1)C 
' " 2 ^tt sin(jC/2) 

Now, we wish to compute the sum d'l and prove that it is negative. For this, 
we need the foUowing trigonometric identity 

cosjk( — cosj(fc — 1)C = (1 — cosjC) cosjk( — sinj/cCsinfcC 

o • ^'^ f ■ ^ • ., > K 

— 2 sm — sm — cos ikC — sm ikC cos — 

2^2 2 
= -2siny sin(A: - i)jC- 
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Therefore, 

n-l 

4:--^sin(fc-l/2)jC. 

j=o 

From the previous lemma, we have that d'l = — cot(^/2), then d'^ is negative 
when ^ — 7r(2fc — l)/n G (0, tt). Therefore, the numbers dj! are negative for 
A; e N n (1/2, (n + l)/2), and so the numbers d'f, decrease for e N n [0, n/2]. 

Since d'g — 2si is positive and d'f. decreases for fc S N n [0, n/2], there can be 
at most one ko G N fl [1, n/2] such that 

4o-i >o>4o- 

Suppose first that fco does not exist. Then d'j, is positive for all k e Nn [0, n/2], 
and dk increases for all e N n [0,n/2]. Therefore, dk > do > for all k e 
Nn [0,n/2], because do — si. 

Now, suppose that fco does exist. Using that Sn-k — Sk, we may prove that 
dk satisfies the equality 

dk-i = (sfc - Sfc-i) = -(s(„_fc)+i - S(„_fe)) -d„^k- (14) 

If n is odd, from the equality ([T^ with k — {n + l)/2, we have that d(^n-i)/2 = 
—d^n~i)/2 =0. If n is even, from the equality (|14p with k = n/2 we have that 
dn/2-i — ^dn/2- Siucc wc did suppose that the fco exists, then 

2(^,1/2-1 = dn/2-l ~ dn/2 = ~'^^l/2 > 0- 

Finally, since d^ increases in Nn(0, fco) and decreases in Nn(fco, [(n— 1)/2]], then 
the numbers dk are positive for fc g {0, [(n — l)/2]}, because dk is positive 
at the end points fc = 0, [(n — l)/2]. In any case, the sums Sk increase for 
fce{0,...,[n/2]}. ■ 
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